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Abstract 

We show that in bounded domains with no-slip boundary conditions, the Navier- 
Stokes pressure can be determined in a such way that it is strictly dominated by 
viscosity. As a consequence, in a general domain we can treat the Navier-Stokes equa- 
tions as a perturbed vector diffusion equation, instead of as a perturbed Stokes system. 
We illustrate the advantages of this view in a number of ways. In particular, we provide 
simple proofs of (i) local-in-time existence and uniqueness of strong solutions for an 
unconstrained formulation of the Navier-Stokes equations, and (ii) the unconditional 
stability and convergence of difference schemes that are implicit only in viscosity and 
explicit in both pressure and convection terms, requiring no solution of stationary 
Stokes systems or inf-sup conditions. 



1 Introduction 

The pressure term has always created problems for understanding the Navier- 
Stokes equations of incompressible flow. Pressure plays a role like a Lagrange 
multiplier to enforce the incompressibility constraint, and this has been a main 
source of difficulties. Our general aim in this paper is to show that the pressure 
can be obtained in a way that leads to considerable simplifications in both 
computation and analysis. 

From the computational point of view, typical difficulties are related to the 
lack of an evolution equation for updating the pressure dynamically and the lack 
of useful boundary conditions for determining the pressure by solving boundary- 
value problems. Existing methods able to handle these difficulties are sophis- 
ticated and lack the robustness and flexibility that would be useful to address 
more complex problems. For example, finite element methods have required 
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carefully arranging approximation spaces for velocity and pressure to satisfy an 
inf-sup compatibility condition jGRj . Projection methods too have typically en- 
countered problems related to low-accuracy approximation of the pressure near 
boundaries |Chl ITe2l IQIDj . Yet much of the scientific and technological signif- 
icance of the Navier-Stokes equations derives from their role in the modeling 
of physical phenomena such as lift, drag, boundary-layer separation and vortex 
shedding, for which the behavior of the pressure near boundaries is of great 
importance. 

Our main results in this article indicate that in bounded domains with no-slip 
boundary conditions, the Navier-Stokes pressure can be determined in a such 
way that it is strictly dominated by viscosity. To explain, let us take Q to be a 
bounded, connected domain in M. N (JV > 2) with C 3 boundary T = dfl. The 
Navier-Stokes equations for incompressible fluid flow in Q with no-slip boundary 
conditions on T take the form 

d t u + W-W + Vp = vAu + / in f2, (1) 
V • u = in n, (2) 
u = on T. (3) 

Here u is the fluid velocity, p the pressure, and v is the kinematic viscosity 
coefficient, assumed to be a fixed positive constant. 

A standard way to determine p is via the Helmholtz-Hodge decomposition. 
We let V denote the Helmholtz projection operator onto divergence-free fields, 
and recall that it is defined as follows. Given any a £ L 2 (fl,M. N ), there is a 
unique q £ H 1 ^) with J Q q — such that Va := a + Wq satisfies 

= / (W)-V0= / (a + X7q) ■ V<£ for all G H^fl). (4) 
Jn Jn 

The pressure p in Q is determined by taking a — u-'Vu — f — vAu. Then Q 
is rewritten as 

d t u + V(u-Vu- f- vAu) = 0. (5) 

In this formulation, solutions formally satisfy c?t(V ■ u) — 0. Consequently 
the zero-divergence condition J5J needs to be imposed only on initial data. Nev- 
ertheless, the pressure is determined from (JjJJ in principle even for velocity fields 
that do not respect the incompressibility constraint. However, the dissipation 
in ||SJ| appears degenerate due to the fact that V annihilates gradients, so the 
analysis of (JSJ is usually restricted to spaces of divergence-free fields. 

Alternatives arc possible in which the pressure is determined differently when 
the velocity field has non-zero divergence. Instead of ©, we propose to consider 

d t u + V(u-Vu- f-vAu) = vV(V ■$). (6) 

Of course there is no difference as long as V • u = 0. But we argue that 
10 enjoys superior stability properties, for two reasons. The first is heuristic. 
The incompressibility constraint is enforced in a more robust way, because the 
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divergence of velocity satisfies a weak form of the diffusion equation with no-flux 
(Neumann) boundary conditions — Due to (J2J , for all appropriate test functions 
4> we have 

/ d t u ■ V4> = v I V(V ■ u) ■ V(/>. (7) 
Jn Jn 

Taking <fi = V • u we get the dissipation identity 

i\ l^-uf + v f |V(V-u)| 2 = 0. (8) 

Due to the Poincare inequality and the fact that J n V • u = 0, the divergence 
of velocity is smoothed and decays exponentially in L 2 norm. Naturally, if 
V • u = initially, this remains true for all later time, and one has a solution of 
the standard Navier-Stokes equations Q-Q- 

The second reason is much deeper. To explain, we recast (0 in the form 
l|T]l while explicitly identifying the separate contributions to the pressure term 
made by the convection and viscosity terms. Using the Hclmholtz projection 
operator V, we introduce the Euler pressure p E and Stokes pressure p s via the 
relations 

V(u-Vu- f) = u-Vu- f + Vp E , (9) 
V(-Au) = -Au + V(V-u) + Vp s . (10) 

This puts JHJ into the form JQ) with p — p E + vp s : 

d t u + u-X7u + Vp E + i/V]3s = vAu + f. (11) 

Identifying the Euler and Stokes pressure terms in this way allows one to focus 
separately on the difficulties peculiar to each. The Euler pressure is nonlinear, 
but of lower order. Since the Hclmholtz projection is orthogonal, naturally the 
Stokes pressure satisfies 

/ |Vp s | 2 < / |Au| 2 ifV-u = 0. (12) 

The key observation is that the Stokes pressure term is actually strictly domi- 
nated by the viscosity term, regardless of the divergence constraint. We regard 
the following theorem as the main achievement of this paper. 

Theorem 1 Let SI c M. N (N > 2) be a connected bounded domain with C 3 
boundary. Then for any e > 0, there exists C > such that for all vector fields 
u G H 2 n Hq (SI, M. ), the Stokes pressure p s determined by \ 1(J\) satisfies 

! \V P s\ 2 <P f \Au\ 2 + C [ |W| 2 , (13) 
Jn Jn Jn 



where (3 = | + e. 
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This theorem allows one to see that (jSJ) is fully dissipative. To begin to 
see why, recall that the Laplace operator A: H 2 (fl) n H^fl) — > L 2 (il) is an 
isomorphism, and note that Vp s is determined by Au via 

Vp s = (I -V - Q)Au, Q — VVA- 1 . (14) 

Equation JSJ can then be written 

d t u + V(u-Vu- f) = v{V + Q)Au 

= vAu-v{I-V-Q)Au. (15) 

Theorem ^ will allow us to regard the last term as a controlled perturbation. 

We can take Am = g arbitrary in L 2 (fl,M. N ) and reinterpret TheorernQJas 
follows. The last term in (|13fl can be interpreted as the squared norm of u 
in i?o(n, K w ), giving the norm of g in the dual space -ff _1 (f2, R"). Thus the 
conclusion of Theorem^is equivalent to the following estimate, which says that 
/ — V is approximated by the bounded operator Q : L 2 (il,M. N ) — » Vi? 1 (fi): 

Corollary 1 For a/Z vector fields g£ L 2 (n,R N ) we have 

\\(I-r-Q)g\\l 2 <P\\g\\ 2 L2 + C\\g\\l^. (16) 

There are several different ways to interpret the Stokes pressure as we have 
defined it. In this vein we make a few further observations. First, note that 
PV(V • u) = for all u in H 2 n R N ), since V ■ u lies in ff^fi). Then 

Vp s = (I-V){Au- V(V-tt)). (17) 

Now „4w := Au — V(V • u) has zero divergence in the sense of distributions and 
is in L 2 (tt, R ), so .Au lies in the space i?(div; 51) consisting of vector fields in 
L 2 (fl,R N ) with divergence in L 2 (f2). By consequence, Ap s = in the sense of 
distributions and so Vp s is in iJ(div; £1) also. By a well-known trace theorem 
(see |GR| . theorem 2.5), the normal components of Au and Vp s belong to the 
Sobolev space i? -1 / 2 (r), and from the definition of V we have 

= / (Vp s - Au) ■ V4> = / <pn ■ (Vp s - .Au) (18) 
Jn Jr 

for all (j> E H^(d). So p s is determined as the zero-mean solution of the Neumann 
boundary- value problem 

Ap s = infi, n ■ Vp s = n ■ (A - VV-)5 on T. (19) 

Furthermore, in two and three dimensions, we have 

Vp s = -(I-V)(\7 x V x u) (20) 

due to the identity V x V x u = — Aw + V(V • u). Green's formula yields 

/ n ■ (V x V x u)<f> = I (V x V x u) ■ V(f> = - / (V x u) ■ {H x V<p) (21) 
Jr Jn Jr 
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and so p s (with zero average) is determined through the weak formulation |JL| 



(Note that V x u G H 1 / 2 ^,^), and nxVcf) £ H-^ 2 {T,R N ) by a standard 
trace theorem GR Theorem 2.11], since V</> lies in if (curl; f2), the space of 
vector helds in L 2 (fl,R N ) with curl in L 2 .) 

As indicated by (|19|1 or 1)22(1 , the Stokes pressure is generated by the tangen- 
tial part of vorticity at the boundary. In the whole space M. N or in the case of 
a periodic box without boundary, the Helmholtz projection is exactly given via 
Fourier transform by V = I — Q and the Stokes pressure vanishes. Essentially, 
the Stokes pressure supplies the correction to this formula induced by the no-slip 
boundary conditions. By consequence, the results of the present paper should 
have nothing to do with the global regularity question for the three-dimensional 
Navier-Stokes equations. But as we have mentioned, many important physi- 
cal phenomena modeled by the Navier-Stokes equations involve boundaries and 
boundary-layer effects, and it is exactly here where the Stokes pressure should 
play a key role. 

The unconstrained formulation JJJJ is not without antecedents in the litera- 
ture. Orszag et al. |OID| used the boundary condition in (|T§)) as a way of enforc- 
ing consistency for a Neumann problem in the context of the projection method. 
After the results of this paper were completed, we found that the formulation JJjJl 
is exactly equivalent to one studied by Grubb and Solonnikov jHSniHS2] . These 
authors also study several other types of boundary conditions, and argue that 
this formulation is parabolic in a nondegenerate sense. They perform an anal- 
ysis based on a general theory of parabolic pseudo-differential initial-boundary 
value problems, and also show that for strong solutions, the divergence satisfies 
a diffusion equation with Neumann boundary conditions. 

Due to our Theorem^ we can treat the Navier-Stokes equations in bounded 
domains simply as a perturbation of the vector diffusion equation dtu = uAu, 
regarding both the pressure and convection terms as dominated by the viscosity 
term. This stands in contrast to the usual approach that regards the Navier- 
Stokes equations as a perturbation of the Stokes system dtu = vAu — Vp, 
V • u = 0. Discussing this usual approach to analysis, Tartar |Ta2l p. 68] 
comments 

"The difficulty comes from the fact that one does not have adequate 
boundary conditions for p. . . .[S] ending the nonlinear term to play 
with /, one considers the Navier-Stokes equations as a perturbation 
of Stokes equation, and this is obviously not a good idea, but no one 
has really found how to do better yet." 

By way of seeking to do better, in this paper we exploit Theorem^in a number 
of ways. In particular, we develop a simple proof of local-in-time existence and 
uniqueness for strong solutions of the unconstrained formulation i|ll|) and conse- 
quently for the original Navier-Stokes equations, based upon demonstrating the 




(22) 



Divorcing pressure from viscosity 



6 



unconditional stability of a simple time-discretization scheme with explicit time- 
stepping for the pressure and nonlinear convection terms and that is implicit 
only in the viscosity term. 

The discretization that we use is related to a class of extremely efficient 
numerical methods for incompressible flow (Til IPel I.TLI IGuS] , Thanks to the 
explicit treatment of the convection and pressure terms, the computation of the 
momentum equation is completely decoupled from the computation of the kine- 
matic pressure Poisson equation used to enforce incompressibility. No stationary 
Stokes solver is necessary to handle implicitly differenced pressure terms. For 
three-dimensional flow in a general domain, the computation of incompressible 
Navier-Stokes dynamics is basically reduced to solving a heat equation and a 
Poisson equation at each time step. This class of methods is very flexible and 
can be used with all kinds of spatial discretization methods | JLj . including fi- 
nite difference, spectral, and finite element methods. The stability properties 
we establish here should be helpful in analyzing these methods. 

Indeed, we will show below that our stability analysis easily adapts to proving 
unconditional stability and convergence for corresponding fully discrete finite- 
element methods with C 1 elements for velocity and C° elements for pressure. It 
is important to note that we impose no inf-sup compatibility condition between 
the finite-element spaces for velocity and pressure. The inf-sup condition (also 
known as the Ladyzhenskaya-Babuska-Brezzi condition) has long been a central 
foundation for finite-element methods for all saddle-point problems including the 
stationary Stokes equation. Its beautiful theory is a masterpiece documented in 
many finite-element books. In the usual approach, the inf-sup condition serves to 
force the approximate solution to stay close to the divergence-free space where 
the Stokes operator "PA is dissipative. However, due to the fully dissipative 
nature of the unconstrained formulation which follows as a consequence 
of Theorem ^ as far as our stability analysis in section HJ is concerned, the 
finite-element spaces for velocity and pressure can be completely unrelated. 

The proof of Theorem ^ will be carried out in section Important in- 
gredients in the proof are: (i) an estimate near the boundary that is related 
to boundedness of the Neumann-to-Dirichlet map for boundary values of har- 
monic functions — this estimate is proved in section |3 see Theorem |2 and (ii) 
a representation formula for the Stokes pressure in terms of a part of velocity 
near and parallel to the boundary. In sectional we also describe the space VS P 
of all possible Stokes pressure gradients (i.e., the range of I — V — Q). In K 3 
it turns out that this is the space of square-integrable vector fields that are 
simultaneously gradients and curls (see Theorem 0] in section Tci.51 below - ) . 

In section 0] we establish the unconditional stability of the time-discretiza- 
tion scheme, and in section El we use this to study existence and uniqueness 
for strong solutions with no-slip boundary conditions. In section El we adapt 
the stability analysis to prove the unconditional stability and convergence of 
corresponding C 1 /C° finite-element methods. 

In section [7| we show that Theorem ^ also allows one to treat the linearized 
equations (an unconstrained version of the Stokes system) easily by analytic 
semigroup theory. We deal with non-homogeneous boundary conditions in sec- 
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tion[S] From these results, in section we deduce an apparently new result for 
the linear Stokes system, namely an isomorphism theorem between the solution 
space and a space of data for non-homogeneous side conditions in which only 
the average flux through the boundary vanishes. 

2 Integrated Neumann-to-Dirichlet estimates in 
tubes 

2.1 Notation 

Let fl C ~M. N be a bounded domain with C 3 boundary T. For any x E fi we let 
&(x) = dist(x, r) denote the distance from x to T. For any s > we denote the 
set of points in Q within distance s from T by 

n s = {xeQ\ < s}, (23) 

and set Q° = f2\f2 s and T s — {x E ft | &(x) = s}. Since T is C 3 and compact, 
there exists so > such that 4> is C 3 in Q So and its gradient is a unit vector, 
with |V$(£)| = 1 for every x E ft So . We let 

H(f) = -V$(x), (24) 

then n(x) is the outward unit normal to T s = d$l c s for s = < E > (a ? ), and n S 

We let (f,g) n = f n fg denote the L? inner product of functions / and g in 
f2, and let || • ||o denote the corresponding norm in L 2 (f2). We drop the subscript 
on the inner product and norm when the domain of integration is understood 
in context. 

2.2 Statement of results 

Our strategy for proving Theorem ^ crucially involves an integrated Neumann- 
to-Dirichlet-type estimate for harmonic functions in the tubular domains tt s for 
small s > 0. 

The theorem below contains two estimates of this type. The first, Q26JI in 
part (i), can be obtained from a standard Neumann-to-Dirichlet estimate for 
functions harmonic in il, of the form 

|(/-m7*)Vp| 2 < \n-W P \ 2 , (25) 

by integrating over r E (0,s), provided one shows that (3q > can be chosen 
independent of r for small r > 0. On the first reading, the reader is encouraged 
to take l|26f) for granted and proceed directly to section 3.2 at this point; it is 
only necessary to replace in the proof of Theorem ^ by the corresponding 

result from (|26|l to establish that the estimate in Theorem ^ is valid for some 
(3 < 1 depending upon Q. 
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The second estimate, in part (ii) , will be used with [3\ close to 1 to establish 
the full result in Theorem ^ for any number (3 greater than |, independent of 
the domain. 

Theorem 2 Let fl be a bounded domain with C 3 boundary, (i) There exists 
(3q > such that for sufficiently small s > 0, whenever p is a harmonic function 
in £l s we have 

A) / {(I-nn^pl 2 < f \n-Vp\ 2 . (26) 

(ii) Let Pi < 1. Then for any sufficiently small s > 0, whenever p is a harmonic 
function in Q s and po is constant on each component ofQ s , we have 

fa f \(I ~ nn^Vpl 2 < f \n-Vp\ 2 + ^ f \p - Po \ 2 ■ (27) 

Our proof is motivated by the case of slab domains with periodic bound- 
ary conditions in the transverse directions. In this case the analysis reduces 
to estimates for Fourier series expansions in the transverse variables. For gen- 
eral domains, the idea is to approximate —A in thin tubular domains Q B by 
the Laplace-Beltrami operator on T x (0, s). This operator has a direct-sum 
structure, and we obtain the integrated Neumann-to-Dirichlet-type estimate by 
separating variables and expanding in series of eigenfunctions of the Laplace- 
Beltrami operator on T. For basic background in Riemannian geometry and the 
Laplace-Beltrami operator we refer to |Au| and | Ta| . 

2.3 Harmonic functions on T x (0, s) 

Geometric preliminaries. We consider the manifold Q = Txl withX = (0, s) 
as a Riemannian submanifold of R N x K with boundary dQ = T x {0,s}. We let 
7 denote the metric on T induced from M. N , let t denote the standard Euclidean 
metric on X, and let g denote the metric on the product space Q. Any vector a 
tangent to Q at z = (y, r) has components ar tangent to T at y and ax tangent 
to I at r. For any two such vectors a and 6, we have 

g{a,b)^j(a r ,b r ) + L(ai,b x ). (28) 

Given a C 1 function z — (y,r) i— » /(y,r) on Q, its gradient Vg/ at z is a 
tangent vector to Q determined from the differential via the metric, through 
requiring 

g(V g f, a)=df -a for all a € T Z Q. (29) 

By keeping r fixed, the function y i— > /(?/, r) determines the gradient vector Vr/ 
tangent to T in similar fashion, and by keeping y fixed, the function r <— > /(y, r) 
determines the gradient vector Vj/ tangent to T. These gradients are also the 
components of Vg/: 



(V e /)r = Vr/, (Vgf)x = V x /« 
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If u = (u 1 , . . . jii"" 1 ) i— ► y = [y , . . . ,y ) is a local coordinate chart for 
r, the metric is given by 7^ du 1 du 3 (summation over repeated indices implied) 
with matrix elements 

Qyk Qyk 

For Icl the identity map serves as coordinate chart. In these coordinates the 
tangent vectors are written (in a form that aids in tracking coordinate changes) 

as . df d df d 

*^ Tjt ^ Qu l du 3 ' Q r q t (^^) 

As usual, the matrix (7^) = (7y )~ 1 - Given two C 1 functions /, / on Q, 

7(V r /,V r /,=7«jV§, l(Vl/ , Vl/) = H. m 

In these coordinates, the (positive) Laplace-Beltrami operators on T and 2 
respectively take the form 

Ar/ = -^(>7«^/), Ax/ = -^/, (32) 



where ^7 = yf det (7^ ) is the change-of- variables factor for integration on T 
if a function / on T is supported in the range of the local coordinate chart then 

f(y)dS(y)= [ f{y{u))^du. (33) 



r 



(Since orthogonal changes of coordinates in R N and K^ -1 leave the integral 
invariant, one can understand ^f 7 ) as the product of the singular values of the 
matrix dy/du.) 

Whenever / e H 1 ^) and / 6 H 2 (T), one has the integration-by-parts 
formula 

J /Ar/ = y\(V r /,Vr/). (34) 

One may extend Ar to be a map from i? 1 (r) — * i/ _1 (r) by using this equation 
as a definition of Ar/ as a functional on iJ 1 (r). In standard fashion |Ta| . one 
finds that 1+ A r : H 1 ^) -> i? _1 ( r ) is an isomorphism, and that (/ + A r ) _1 is 
a compact self- adjoint operator on L 2 (T), hence L 2 (T) admits an orthonormal 
basis of eigenfunctions of Ar- Since the coefficient functions in i.Y2li are C 1 , 
standard interior elliptic regularity results C |GTI Theorem 8.8], |Tal p. 306, 
Proposition 1.6]) imply that the eigenfunctions belong to H 2 (T). We denote 
the eigenvalues of Ap by v 2 , k = 1,2, . . ., with = v\ < v% < . . . where Vf, — > 00 
as k — > 00, and let tpk be corresponding eigenfunctions forming an orthonormal 
basis of L 2 {Y). If Ar^ = then ip is constant on each component of T, so if rn 
is the number of components of T, then = v m < v m +i- 

In the coordinates u = (u,r) 1— > z = (y,r) for Q, the metric g takes the 
form jij du 1 du 3 + dr 2 , and the Laplace-Beltrami operator Ag = Ar + Ax. 
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Similar considerations as above apply to Ag, except Q has boundary. Whenever 
/ € Hl{Q) and / £ H 2 (G) we have 

f /Ag/= [ g(Vgf,Vgf). (35) 
Jg Jg 

One extends Ag to map H 1 (G) to H^ 1 (G) by using this equation as a definition 
of Agf as a functional on Hq(Q). 

We introduce notation for L 2 inner products and norms on Q as follows: 

(f,f)g=[ff ll/lle = / l/l 2 , (36) 

Jg Jg 

(Vr/,V r /)g = / 7 (Vr/,Vr/), ||V r /|| g = / 7(Vr/,V r /), (37) 
Jg Jg 

(Vl/,Vl/) S = / (d r f)(d r f), ||Vx/||g= / (9 r /) 2 , (38) 

Jg Jg 

(Vg/, Vg/)g = / , 9 (Vg/, Vg/) = (V r /, V r /)g + (Vj/, Vx/) S) (39) 
Jg 

l!Vg/||g = / 5 (Vg/, Vg/) = ||V r /||g + ||Vl/|||. (40) 

Jg 

Lemma 1 Suppose f £ H 1 {G) and Agf = on C*=rx2" where X — (0, s). 
Then, (i) there exists (3q S (0, 1) independent of f such that 

A>||V r /||g < llVi/Hg, (41) 

and (m) 

l|V r /||g<||Vx/||g + ^||/-/o||g (42) 

whenever fo is constant on I\ x (0, s) /or every component Ti ofT. 

Proof: Suppose Agf = on 5. Since the coefficient functions in (|32|) are 
C 1 , the aforementioned interior elliptic regularity results imply that that / £ 
H? oc (G)- For any r £ (0, s), fixing r yields a trace of / in i? 1 (r), and as a function 
of r, we can regard / = /(y, r) as in the space L 2 ([a, 6], H 2 (T))(lH 2 ([a, b],L 2 (T)) 
for any closed interval [a, 6] C (0,s). Now, for each r we have the L 2 (r)- 
convergent expansion 

f(y,r) = J2Hk,r)My) (43) 

k 

where 

f(k,r) = J f(y,r)Mv)dS(y). (44) 
For each k £ N, the map r i— > /(fc, r) is in H 2 OC (0, s) and 

dj(k,r) = J d r f(y,r)i> k (y)dS(y). (45) 
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For any smooth £ £ Cq°(0, s), taking f(y,r) — tp k (y)£,(r) we compute that 

Vr/ - e(r)V r ^ fc , d r f = (46) 
and so by (|3*5j) . (|2"5j) . and (|3~4*|) . we have 

= J(A s f)f = JJ^ ( 7 (V r /, V r /) + (d r f)(d r fj) 
f(r) / 7(Vr/,VrVfc)+ / (ft-fl / (ft./)Vfc 



= / £(r) y /A r Vk + J^d r Od r f(k,r) 

(t{r)4f{k, r) + (d r Od r f(k, r)) dr. (47) 

Therefore /(fc, •) is a weak solution of d 2 / = z^/ in i?j^ c (0, s) and hence is C 2 
and it follows that whenever v k ^ 0, there exist a k , b k such that 

/(fc, r) = a k sinh z^r + b k cosh i/^r, T — r — s/2. (48) 

Now 

ll/lll=E [ S \f(k,r)\ 2 dr, (49) 
fc • / ° 

H v r/lle = E [ S WJ(k,r)\ 2 dr, (50) 

HVi/||| = E f S \drf(k,r)\ 2 dr. (51) 
fc J ° 

Let 7^ = smn2 v k r dr. Then 7^. increases with fc, and 

/•s/2 rs/2 U 2 S 3 

7fe + s= / cosh 2 v k Tdr> I (l + ^ 2 r 2 )dr> . (52) 
.7-8/2 J-s/2 12 

Whenever i/fc ^Owe get 

\f(k,r)\ 2 dr = \a k \ 2 lk + \b k \ 2 ( lk + s), (53) 

|9 r /(fc, r)| 2 dr = ^ 2 (|a fc | 2 (7fc + «) + M 2 7fc), (54) 
and since /3o(7fc + s ) < 7fc where /3q = 7m+i/( 7 m+i + s), it follows 

A) P Wkf(k,r)\ 2 dr< f \d r f(k,r)\ 2 dr, (55) 
Jo Jo 

f W k f(k,r)\ 2 dr< f \d r f(k,r)\ 2 dr + ^ f \f(k,r)\ 2 dr. (56) 
Jo Jo s Jo 

The results in (i) and (ii) follow by summing over k. □ 
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2.4 Global coordinates on T x (0, s) 

It will be important for comparison with the Laplacian on fl s to coordinatize Q 
for small s > globally via the coordinate chart Cl s — » Q given by 

x h-> z = (y, r) = (a; + $(x)n(a;), (57) 

In these coordinates, the metric on Q that is inherited from R^" 1 " 1 has the 
representation <?y dx l dx 3 with matrix elements given by 

dz k dz k dy k dy k dr dr 
^ dx 1 dx 3 dx 1 dx 3 dx % dx 3 



Let us write di = d/dx l and let V/ = {d\f,..., On/) denote the usual gradient 

%rij = dji 



vector in R . The components of ft are rii = —9;$ and so diUj = 9,7?,;, meaning 



the matrix Vn is symmetric. Since |n| 2 = 1 we have riidjni = = n^rij. Then 
the N x N matrix 

^ = I - nrt* + $Vn = (J — rari*)(J + $Vrt)(I - rm*), (59) 
ax 

and the matrix 

G = Oij) = (I- nn l ){I + $Vn) 2 (I - nn l ) + tin 1 = (I + $Vn) 2 . (60) 



With y/g = y/dct G, the integral of a function / on Q in terms of these 
coordinates is given by 

r f= I fVddx. (61) 



Given two G 1 functions /, / on Q, we claim that the following formulae are 
valid in the coordinates from (|5T|) : 

3(V e /, Vgf) = (VfYG-^Vf) = g^djdjf, (62) 
7(Vr/,Vr/) - (V/)*(J - fHP)G-\l - nn*)(V/), (63) 
t(Vi/,Vi/) = (n-V/)(n-V/) = (V/)*rfn*(V/). (64) 

Of course (|62|l simply expresses the metric in the ^-coordinates from 157|) . To 
prove H64f> . first note that along any curve r i— > a;(r) satisfying 9 r 2; = n(x) we 
find d T n(x) — rijdjrii = 0, so n(cc) is constant and the curve is a straight line 
segment. Hence in the chart from (|57ll . n(x) = n(y) and we have x = y — rn(y). 
Given a G 1 function / then, we find that in these ^-coordinates, 

d r f{y,r) = (d r Xj)(djf) = n^djf = n-V/, (65) 

and (Hll) follows from (J2U- Finally, (HSJ) follows directly from (JHIJ) and (M)) 

using l)28fl — since nn'Vn = we have nn*G = nn* so nn* = nn*G _1 and 
hence 

(7 - r?n*)G _1 (/ - rm*) = G" 1 - nrt*. (66) 
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2.5 Proof of Theorem H 

Let (3\ < 1. Suppose Ap = in fi s . We may assume p € -ff 1 (fi s ) without loss 
of generality by establishing the result in subdomains where $(a;) € (a, 6) with 
[a, 6] C (0, s) and taking a — > 0, b — > s. We write 

P = Pl +P2, 

where pi G Hq(CI s ) is found by solving a weak form of Agpi — Agp: 

(Vffpi, V e 0) e = (V e p, V e 0) g for all e (67) 

For small s > 0, G = (g lj ) = I + 0(s) and ^/g = 1 + O(s). Since (Vp, Vpi) = 0, 
taking <f> = p\ we have 

l|V s pi|||= / (Vp)*(G- 1 V5-/)Vpidic<Cs||Vp||n.||V ( ;pi||t; 1 (68) 

where C is a constant independent of s. By Poincare's inequality we also have 

\\pi\\l<^gPi\\ 2 g (69) 



7T 

since the eigenvalues of Ag on the product space T x [0, s] with Dirichlet bound- 
ary conditions all have the form /i = ^+j 2 7r 2 /s 2 for j, k £ N, so that /x > tt 2 /s 2 . 

Let us first prove part (ii). For < e < 1, using i jSSjl . IpxT]) and i|3T|) we 
deduce 

j|(/-rH<)Vp|| 2 s < (l + Ca)||Vrp||| 

< (1 + Cs) ((1 + E )||VrP2||| + (1 + e-^llVrPill^) 

< (1 + Cs)(l + e) (HVrpallS + e^CVU Vp|| 2 J . (70) 

Now P2 — P ~ Pi satisfies Agp 2 — in f2,j and P2 € H 1 (Q), hence for any p 
constant on each component of f2 a we have 

HVrPalle < \\Vi P 2\\ 2 g + ^\\p2-po\\ 2 g, ( 71 ) 
||ViP2||| < (l+e)||Vip|| 2 +(l+e- 1 )||Vip 1 || 2 

< (l + s)(l + Cs)(\\n-Vp\\ 2 ns +e- x C 2 s 2 \\Vp\\ 2 na ), (72) 

^||pa-Po|li < §2 (lb-Po||S + IIPilll) 



< ^l|p-po||| + ^||V e pi||| 
24 

< -(l + Cs)\\p~ Pa \\l s +C 2 s 2 \\Vp\\ 2 ns (73) 
Presuming Cs < -|e, assembling these estimates yields 

|| (7 -rm^Vp^ < (l + £ ) 4 ( ||n.Vp|| 2 a + -Hp-poll^ +£ || V p| |n 



< (1+e 



) 5 (lln-Vpll^ + ^\\p~Po\\l s + e\\(I-Hn t )^pf n ^ , (74) 
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since |Vp| 2 = \n-Vp\ 2 + \(I—nn t )Vp\ 2 . Fixings > small so that (l + e)~ 5 — e > 
Pi proves part (ii). 

To prove part (i), instead of l(7T|) we use 

/3o||V r p 2 ||2 < HVzPall! (75) 
(from part (i) of Lemma ^1 together with (|70|l and i|72|) and obtain 

hW-m^vpWk < (i + £) 4 (ll^-v P ||^ + |||v P ||^) 

< (l+ef(\\n-Vp\\l s +e\\{I-m! t )Vp\\l s ). (76) 
Now taking e > so small that e(l + e) 5 < /3o finishes the proof. □ 

3 Analysis of the Stokes pressure 

The main purpose of this section is to prove Theorem ^ We also describe the 
range of the map u i— > Vp s from velocity fields to Stokes pressure gradients. 
For motivation for the proof of Theorem ^ the reader can proceed directly to 
section T3. 31 at this point. Here, we first establish some key preliminary results. 

3.1 An L 2 estimate 

The following L 2 estimate on the Stokes pressure will be used to obtain the full 
result of Theorem ^ for arbitrary (3 > | . It is not needed to prove the weaker 
statement that (I13[) holds for some j3 < 1. 

Lemma 2 Let fl C M. N (N > 2) be any bounded connected domain with C ' 
boundary. For any s > 0, there is a constant C > so that for any u € 
H 2 nH^(Q,R N ), the associated Stokes pressure p s defined by with zero 
average satisfies 

||p fl ||<e||AiZ||+C||tZ||. (77) 
Proof: For any <j> £ L 2 (fl), define ip by 

Aip = </> - n-V?A| r = 0, (78) 

where 4> is the average value of <f> over f2. Recall p s — 0. Then, 

(p s ,(/>) = (p s , <j> -4>) = ( Ps ,AtJj) = -<Vp s , W>). (79) 

From J22J), we know (Vp s ,Vip) = (V x u,n x VV>) r when N — 2 or 3. For 
general TV, using the notation di :— d/dxi and automatic summation upon 
repeated indices, from (|19|l we derive 

(Vp s , VV>) = i J (djUi - d i Uj){n j d i 'tl) - n.djtf;). (80) 
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Plug (|8U[) into H79[). take the absolute value and use the trace theorem to get 

\(ps,4>)\ < coll Vw|| jL 2 ( r) || V-0|| i2(r) < c 1 \\u\\ H 3/2 m \\ip\\ H3/ 2 in) . (81) 
By the regularity theory for Poisson's equation lf?%|) , 

U\\ H v H n) <c 2 H- 0|| <caW. (82) 
By a standard interpolation theorem, for any S > 0, there is a constant c, so 

N|j*»/a ( n) <*||AtZ||+c||tf||. (83) 
Plugging and into i|5T)l. we get 

|<Ps,0)| < (<5||At?||+c|H|) Cl c 2 ||0||. (84) 

Thus, 

\\ Ps \\ = sup %^ < e||Au|| + cc lC2 \\u\\. □ (85) 

11011 

3.2 Identities at the boundary 

A key part of the proof of Theorem ^ involves boundary values of two quantities 
that involve the decomposition of u = (I — nn t )u + nrt'u into parts parallel and 
normal to the boundary. Our goal in this subsection is to prove the following. 

Lemma 3 Let ft C I 1 *' be a bounded domain with boundary T of class C 3 . 
Then for any u € H 2 (Vl,M. N ) with u\r = 0, the following is valid on T: 

(i) V • ((/- =0 in H l / 2 {T). 

(ii) ft - (A — VV-)(nn*w) =0 in H- 1 / 2 ^). 

The proof will reduce to the case u <E C 2 (Cl,M. N ), due to the following density 
result. 

Lemma 4 Let f2 C R N be a bounded domain with boundary T of class C 2 ' a 
where < a < 1. Then for any u £ H 2 n Hq(£1,M. n ) , there exists a sequence 
u*k G C 2,a (Cl) such that u k \r = and \\u k — ^\\h 2 (Q) — * 0. 

Proof: Define / = Ait. Since / € L 2 , we can find a sequence f k £ C 1 (f2) so 
that ||/fc — /|| jji — > 0. Construct by solving 

Au k = fk, u k \ gn = 0. 

Classical elliptic regularity theory in Holder spaces (see jCTj . theorem 15.13) 
says that a unique u k exists and is in C 2,a (f2). By standard regularity theory 
in Sobolev spaces, 

\\u k - u\\ H * < C\\f k - f\\ L * - 0. □ 
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Proof of Lemma |3J To begin, recall ft — — V<i>. Equality of mixed partial 
derivatives yields djUi — diUj for all i,j — 1,...,N. Together with the fact 
riiTii = 1, we infer that for small s > 0, throughout £l s we have 

riidjUi = and nidiUj = 0. (86) 

(i) First, for any / G C (CI), if / = on T then V/ || ft on T, which means 

(I - ?m*)V/ = 0, or {d k - nunfij)] = for fc = 1, . . . , N. (87) 

Now suppose u G C 2 (f2, R N ) with u = on T. Then, after taking derivatives in 
fl s for some s > and then taking the trace on T, using (|87() we get 

V • ((/ — nn')u) = dj (Uj — rtjUkUk) = djUj — rijUkdjUk = 9jUj — dkUk = 0. 

For general u 6 H 2 (Sl,R N ) with u\r = 0, the expression V • ((/ — nn^u) is 
in i? 1 (fi s ) for small s > and hence is in if 1 / 2 (r) by a trace theorem. After 
approximating u using Lemma 0] we obtain the result in (i). 

(ii) At first we suppose u £ C 2 (f2,R JV ) with u = on T. We claim in fact 
that for any / G C 2 (Cl) with /| r = 0, 

n- (A- VV-)(n/) = onT. (88) 

This yields (ii) by taking / = ft - u. We prove l|88|l in two steps. 

1. The formula in (i) holds in C(T) if u is C 1 . Since / — rin 1 = (I — riri*) 2 , 
we can use u = (I — nn*)V/ in (i) to find that 

V • ((J - nn^Vf) =0 on T. (89) 

2. Using (|86J) it is easy to verify the following identities in fi s : 

n-A{nf) = Af + fn-An, (90) 
n-VV-(n/) = (m 1 ) : V 2 / + (V ■ n)n • V/ + /n ■ VV • n, (91) 
V-(nn*V/) = (nn ( ) : V 2 / + (V ■ n)n • V/. (92) 

Here (rm*) : V 2 / := niUjdidjf. It directly follows that 

ri- (A - W-)(nf) = V • (J - wj*)V/ + /n • (A - VV-)«. (93) 

Using this with l|89|l proves Ij88(l . and establishes (ii) when u G C 2 (f2) with u = 

on r. 

To establish (ii) for arbitrary u G ff 2 (J7,R A '), we restrict to f2 s for small 
s and let a = (A - VV-)(wH t M)- Then a G L 2 (fi s ,M JV ) and V ■ a = in the 
sense of distributions, so a G i?(div;f2 s ) and a well-known trace theorem (see 
|GR| . theorem 2.5) yields that the map H 2 (Vt s ,M. N ) -> iJ(div;O s ) iJ- 1 / 2 (r) 
given by u i— * a i— * n • a is continuous. To conclude the proof, simply apply the 
approximation lemma above to infer ft ■ a|r = 0. □ 
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3.3 Identities for the Stokes pressure 

Given u £ 77 2 D H$(Cl,R N ), recall that 7>(VV • u) = 0, so that the Stokes 
pressure defined in (jlOJI satisfies 

Vp s = Au - VV • u - 7>Au = (7 - V){A - W-)u. (94) 

Also recall that whenever a £ L 2 {n,R N ) and V • a £ L 2 (Q), n- a £ 77 _1 / 2 (r) 
by the trace theorem for H(div; fl). If V • a = and n ■ a\r = 0, then we have 

V(^>) = for all (f> £ Tf^fi) and this means (I - V)a = 0. Thus, the Stokes 
pressure is not affected by any part of the velocity field that contributes nothing 
to n ■ a\r where a = (A — VV-)u. Indeed, this means that the Stokes pressure 
is not affected by the part of the velocity field in the interior of f2 away from 
the boundary, nor is it affected by the normal component of velocity near the 
boundary, since n ■ (A - VV-)(nn*w)| r = by Lemma El 

This motivates us to focus on the part of velocity near and parallel to the 
boundary. We make the following decomposition. Let p : [0, oo) — > [0, 1] be a 
smooth decreasing function with p(t) = 1 for t < | and pit) = for t > 1. For 
small s > 0, the cutoff function given by £(x) = p($>(x)/s) is C 3 , with £ = 1 
when $(a;) < and £ = when <i>(x) > s. Then we can write 

u = u± + M|| (95) 

where 

u± = (1 — £)tt + £nn*u, = £(7 - nn t )u. (96) 

Since uj_ = (rm 4 )u in O s /2, with aj_ = (A — VV')mi we have 

nieL 2 ^!"), V-a ± =0 and n-aj r =0 (97) 

by Lemma|3Jii). Hence (d±,'V(f)) = for all <f> £ 77 1 (11) , that is, 

(7-P)(A- VV-)#l = 0. (98) 

Combining this with (|94|l and (|95|) proves part (i) of the following. 

Lemma 5 Let Vl C R 6e a bounded domain with C 3 boundary, and let u £ 
H 2 n H^iVl, M. N ). Let p s and 6e defined as in \94\ ) and Mb]) respectively. 
Then 

(i) The Stokes pressure is determined by u\\ according to the formula 

Vp s =(7-P)(A-VV-)w||. (99) 

(ii) For any q £ 77 1 (17) that satisfies Aq = in the sense of distributions, 

(At?||-Vp flJ Vg> = 0. (100) 
(Hi) In particular we can let q — p s in (ii), so (Au\\ — Vp s , Vps) = and 

||AiT|||| 2 = ||Au|| - Vp.,|| 2 + ||Vp.,|| 2 . (101) 
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Proof: We already proved (i). For (ii), note by Lemma |3fi) we have 

V • t*|| | r = 0, (102) 
so V ■ G flo(fi), thus (VV ■ U||, Vq) = -<V ■ u\\,Aq) = 0. Now (i) entails 

(Vp s ,Vg> = (A5||,Vg>. (103) 
This proves (ii), and then (iii) follows by the L 2 orthogonality. □ 

3.4 Proof of Theorem H] 

Let e > and (3 = § + e. We fix ft < 1 such that 1 + e := ftl + > 1, 
and fix s > small so Theorem (ii) applies in fl s with this ft . Let u € 
-ff 2 n -ffp^jR^) and define the Stokes pressure Vp s by (|10|l and the decompo- 
sition u = u± + uii as in the previous subsection. Then by part (iii) of Lemma[S] 
we have 

||Au|| 2 = ||Au_ L || 2 + 2(At? ± ,Au||) + ||Au|| - Vp s || 2 + || Vp s || 2 . (104) 

We will establish the Theorem with the help of two further estimates. 

Claim 1: For any e\ > 0, there exists a constant C\ > independent of u such 

that 

(Aui,Au||) > -£i||Au|| 2 - Ci||Vu|| 2 . (105) 
Claim 2: For any e\ > there exists a constant C 2 independent of u such that 

||AtT|| - Vp s || 2 > ^||Vps|| 2 - ei || Auf - C 2 \\Vu\\ 2 . (106) 

Proof of claim 1: From the definitions in H96fl . we have 

Au± = £nn t Au + (1 - £)Au + R u M\\ = f(I - nH^Au + R 2 , (107) 

where ||i?i|| + \\R2W < C||Vu|| with C independent of u. Since / — nn* = 
(I -fin*) 2 , 

(£nn*Au + (1 - £)Atf) ■ (f(J - nn*)Att) = + f(l - - nri*)A<u| 2 > 0. 

This means the leading term of ^Attj_, Ait||) is non-negative. Using the in- 
equality \(a,b)\ < (ei/C)||a|| 2 + (4C/ei)||6|| 2 and the bounds on i?i and R 2 to 
estimate the remaining terms, it is easy to obtain I|1U5|) . 
Proof of claim 2: Recall that u\\ is supported in f2 s , and note 

Au|| = £(J - nn 1 ) Au + R 3 (108) 

where H-R3II < C||Vu||. Since n ■ (I — nn t )Au = we find 

\\ft- Affylln. <C 2 ||Vtt|| (109) 
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with C 2 > independent of u. We use \a + b\ 2 > (1 — £^2 ) I £*| 2 — l a | 2 / £ 2 to get 

||Au||-V Ps ||^ > / |V Ps | 2 + f \n- (Au|, - Vp s )| 2 

> / |V Ps | 2 + (l-£ 2 ) / \n-Vp s \ 2 --( |n-AS||| 2 . (110) 

Next we use part (ii) of Theorem |21 with po = and with /3i J n |n • Vp| 2 added 
to both sides, together with Lemma|21and Poincare's inequality, to deduce that 

V / |VPs| 2 </ \n-Vp s \ 2 +eJ \Au\ 2 + c[ |W| 2 . (Ill) 

Taking 1 — e 2 = A and combining I1U9I) . 111UI) and l|lll|) establishes Claim 2. 

Now we conclude the proof of the theorem. Combining the two claims with 
(Hull , we get 

(1 + SeOIIA^II 2 > (l + y) IIVPsII 2 - (C 2 + 2C 1 )\\Vu\\ 2 . (112) 
Multiplying by (3 and taking e% > so that 3ei < Eq concludes the proof. □ 
3.5 The space of Stokes pressures 

According to (|21|I - H22|I . the space of Stokes pressures, obtainable via l|10[l from 
velocity fields u E H 2 n Hq(Q.,M. n ), can be characterized as the space 

S p := {pe H\Q)/R | Ap = in fl and n-Vp| r G Sr}, (113) 

where Sr is the subspace of ff _1 / 2 (r) given by 

S r := {/ = n • (A - VV-)u|r I u G H 2 n H^{Q,R N )}. (114) 

The Stokes pressure p with zero average is determined uniquely by f — n-Vp|r G 
Sr, with ||p||ifi(n) < C||/llif- 1 /2(r) by the Lax-Milgram lemma. 
The space Sr may be characterized as follows. 

Theorem 3 Assume O C M. N is a bounded, connected domain and its boundary 
r is of class C 3 . Denote the connected components of T by Ti, i = 1, . . . ,m. 
Then 

S r = {fe H-^ 2 (T) | J f = fori = l,...,m}, 

and moreover, the map u i— ► ft • (A — VV-)t*|r /rom iJ 2 n -ffg (^i K ) to Sr 
admits a bounded right inverse. 

Proof. First we check the necessity of the integral conditions. Let u <G 
H 2 n H^(Q,,R N ) and let / = ft ■ (A - VV-)u| r . For each connected compo- 
nent Ti of r, there is an Sj > small enough and a smooth cut-off function pi 
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defined in £1 which satisfies Pi(x) — 1 when dist(x,I\) < Sj and pt(x) — when 
dist(a;,r i ) < s l for all j ^ i. Let a = (A - VV-)(piU). Then a € L^fi,^) 
and V ■ a = 0, so 

/ /= / H-a= I V-a = 0. (115) 

Next, let / e i7 _1 / 2 (r) with J r / = for all i. Treating each boundary 
component separately, we can then solve the problem 

A r ip = -f onT, J ^ = fori = l,...,m, (116) 

where Ap is the (positive) Laplace-Beltrami operator on T. Denote the mapping 
/ <->■ ip by T. Then T: H- 1 (Y) ~> ^H 1 ") is bounded (|Jul theorem 1.71, theo- 
rem 4.7], [H| p. 306, Proposition 1.6]). Also T: L 2 (T) -> £T 2 (r) is bounded, by 
elliptic regularity theory |Tal p. 306, Proposition 1.6]. So, interpolation implies 
(see [EMI vol I, p. 37, Remark 7.6]) 

HV'IU^(r) < C\\f\\ H - 1/HT y (117) 

Now by an inverse trace theorem |RRI Theorem 6.109], there exists a map 
ip h-> g g if 3 (ft) with 

q = and n-Vq^ip on T, ||?||i?3(n) < C||^|| ff 3/2 (r) . (118) 

We may assume q is supported in a small neighborhood of T. Define 

u={I-nn t )Vq. (119) 

Then / i-> u is bounded from 5 r to iJ 2 n ^(Q.M^). We claim 

ft - (A — W-)u = f onT. (120) 

The proof of this claim amounts to showing, by calculations similar to those 
in the proof of Lemma[31 that the normal derivative n-V and normal projection 
nW commute on the boundary with the tangential gradient and divergence 
operators (I — wj*)V and V • (J — nn*) for the functions involved. 

Pirst, since ft ■ u — 0, by expanding A(n ■ u) we get 

ft- Au = -(An) ■ u- 2Vn : Vu = onT, (121) 

since for each i, Vrii is tangential and S/ui is normal to T — indeed, using 
djUi = diUj and (|86|l and (|87|l . we have that 

Vn : Vm = (djTii){djUi) = (dinj)(rijnkdkUi) = on T. (122) 

Next we calculate in ft that 



n ■ VV ■ u = V • (n • Vu) - Vn : V«. 



(123) 
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Note that ft ■V(ftft t ) — by JSfJJ), so ft ■ V commutes with I — ftft 1 in f2. Then 
since u = (I — ftft^u from (|119fl we get 

ft-Vu = (I — nn f ')(n • V)u = (J — nn t )(n • V)Vq. (124) 

Now 

(n-V)Vq = V(n-Vg) - a (125) 

where 

a» = {d i nj){d j q) = (djrii)(djq) (126) 

This quantity lies in iJ 2 (fl) and vanishes on F since Vq = (nn*)Vg on L. (This 
can be proved by approximation using Lemma ^) Using part (i) of Lemma [31 
we have that V • (/ — rm*)a = on L. Combining 11211) (|125f) we conclude that 

ft- (A- VV-)S= -V- (I -fin^Vin-Vq) on L. (127) 

But it is well known that at any point x where $(x) = r£ (0, s), for any smooth 
function <f) on il s , 

V • (7 - rm*)V<£ = A<£- (V ■ n)(n-V<£) - (n-V) 2 </> = -Arv(</>|r r ). (128) 

where Ar r is the Laplace-Beltrami operator on L r . So taking r^Owe see that 
the right hand side of (|127[) is exactly — A r (n-Vg|r)- So by l|lltj|) and !|118|) we 
have established the claim in i|120[l . This finishes the proof. □ 

Remark 1. Given a velocity field u € H 2 n ^(OjM 3 ), the associated Stokes 
pressure is determined by the normal component at the boundary of the curl 
of the vorticity uj — V x u, which is a vector field in 7J 1 (r2,R 3 ). A question 
related to Theorem|2]is whether the space Sr of such boundary values ft ■ V x ui 
is constrained in any way, as compared to the space of boundary values ft ■ V x v 
where v € 7J 1 (r2,R 3 ) is arbitrary. 

The answer is no. In [Tell Appendix I, Proposition 1.3], Temam proves 

V x iJ^^R 3 ) = {g e L 2 {n,M. 3 ) | V ■ g = 0, J ft-g = 0\/i}. (129) 

Clearly S r C ft ■ V x ii^^R 3 ) by ifTHjl . For the other direction, let v € 
F 1 (0,IR 3 ) be arbitrary, and let / = ft ■ V X v\r- By Ijl29jl or otherwise, / € 
i?- 1 /2(r) and J r . / = for all i, hence / G S r . This shows that for N = 3, 

Sr =r?< V x iJ^R 3 ). (130) 

A related point is that for TV = 3, the space of Stokes pressure gradients 
V<S P can be characterized as the space of simultaneous gradients and curls. 

Theorem 4 Assume ft C R 3 is a bounded, connected domain and its boundary 
r is of class C 3 . Then 

vs p = VH\n) n v x /^(^r 3 ). (131) 



Divorcing pressure from viscosity 



22 



Proof. Indeed, V5 P C V x H 1 by (Tg^ and Theorem^ On the other hand, if 
g = V xv = Vp then Ap = V • g = and ra -Vp|r G Sr by |JT25|| and Theorem^ 
so Vp G VS P . □ 

Remark 2. In the book |Telj (see Theorem 1.5) Temam establishes the or- 
thogonal decomposition L 2 (fl,M. N ) = H © H\ © H2, which means that for any 
g€L 2 (Q,R N ), 

g = Vg + Vq + VA^V-s, (132) 

where g satisfies Aq = and n ■ Vg|r — n ■ (g ~ VA _1 V • g). By contrast, we 
have shown 

g = Vg + Vp + VV • A" 1 *? (133) 

where p satisfies Ap = and n ■ Vp\r = n ■ (g — VV ■ A -1 <?), i.e., p is the Stokes 
pressure associated with A -1 *?. Thus the map g 1— > Vp — Vg is the commutator 
VA _1 V-— VV- A -1 . The decomposition (|132|l is orthogonal, and q satisfies in- 
Vq, l) r = 0. In our decomposition l|133(l . the gradient terms are not orthogonal, 
but the Stokes pressure term enjoys the bounds stated in Corollary^ an d if T 
is not connected, it has the extra property that (n ■ Vp, l) r = for every i. 

4 Unconditional stability of time discretization 
with pressure explicit 

In this section we exploit Theorem ^ to establish the unconditional stability of 
a simple time discretization scheme for the initial-boundary-value problem for 
JSJ, our unconstrained formulation of the Navier-Stokes equations. We focus 
here on the case of two and three dimensions. In subsequent sections we shall 
proceed to prove an existence and uniqueness theorem based on this stability 
result. 

Let f2 be a bounded domain in M. N with boundary T of class C 3 . We consider 
the initial-boundary-value problem 

dtti + u-Vu + Vpv + vVps =vAu + f (t>0, aiefi), (134) 

u = {t > 0, x eT), (135) 
u = u in (t = 0, x en). (136) 

We assume u in G H^{Q, R N ) and / G L 2 (0, T; L 2 {Q, R N )) for some given T > 0. 
As before, the Euler and Stokes pressures p E and p s are defined by the relations 

V(u-\7u- f) = u-Vm - / + Vp E , (137) 
V{-Au) = -Au + V(V • u) + Vp s . (138) 

Theorem ^ tells us that the Stokes pressure can be strictly controlled by 
the viscosity term. This allows us to treat the pressure term explicitly, so 
that the update of pressure is decoupled from that of velocity. This can make 
corresponding fully discrete numerical schemes very efficient (see |.TL| ) . Here, 
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through Theorem ^ we will prove that the following spatially continuous time 
discretization scheme has surprisingly good stability properties: 

vAu rl+1 = f n - vJ 1 ■ Vu n - Vpl - ^Vp™, (139) 



At 

Vp% = (I-V)(f n -u n -Vu n ), (140) 

Vp™ = (I-V)Au n - V(V-tt™), (141) 

u"| r = 0. (142) 



We set 



1 



(n+l)A* 



At inAt 

and take it e # 2 n H£(£l,R N ) to approximate M in in H^(Q,M N ). It is evident 
that for all n = 0, 1, 2, . . ., given tf e if 2 n Hq one can determine Vp^ <= ^ 2 
and Vp™ G L 2 from (|140|l and (|141|l and advance to time step n + 1 by solving 
<|139|) as an elliptic boundary-value problem with Dirchlet boundary values to 
obtain u n+1 . 

This simple scheme is related to one studied by Timmermans et al. |Tij . In 
the time-differencing scheme described in |Ti| for the linear Stokes equation, the 
pressure p n = p™ + vp™ is updated in nearly equivalent fashion, if one omits the 
velocity correction step that imposes zero divergence, and uses first-order time 
differences in (15) and (18) of |Tij. Also see |Fel IHIiSl UT\ . 

Let us begin making estimates — our main result is stated as Theorem 
below. Dot UnS} with -Au n+1 and use (fl4T1|l and \\I - V\\ < 1 to obtain 



— (II 
2AtV 



w" +1 || 2 - ||v?r|| 2 + ||v?r +1 - Vu"|| 2 ) + v\\Au n+1 2 

< ||Au" +1 || (2\\f n - u 11 • W"|| + H|Vp™ I 



< %||Au" +1 || 2 + — \\f n - • Vu"|| 2 + J(||Au" +1 || 2 + ||Vp™|| 2 ) 
2 E\ 2 

(144) 

for any E\ > 0. (This is not optimal for Vp™ but is convenient.) This gives 
i-(|| W" +1 || 2 - || W"|| 2 ) +(u- ei )|| A^f 



< f (Jl/T + K l • ViT|rJ + H|Vp™ II- (145) 
Fix any j3 with | < j3 < 1. By Theorem Q] one has 

H|Vp™ || 2 < AiT|| 2 + vCpfJiPf. (146) 
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Using this in ^ 145(1 . one obtains 

i-(||W" +1 || 2 -||Vtr|| 2 ) e 1 )(||AtT+ 1 || 2 - IIA^H 2 ) 

+ (//-ci-^/3)||A^|| 2 

< — (ll/"l| 2 + IK • Vu n || 2 ) + iyC( 3 \\Vu' l \\ 2 . (147) 

At this point there are no remaining difficulties with controlling the pres- 
sure. It remains only to use the viscosity to control the nonlinear term. We 
focus on the physically most interesting cases N = 2 and 3. We make use of 
Ladyzhenskaya's inequalities |La| 

? 4 <2(7 Ail |Vg| 2 ) (N — 2), (148) 

? 4 <4(7 9 2 ) 1/2 ([ |V 5 | 2 V A " (N = 3), (149) 

valid for g £ iJ 1 (R Ar ) with N = 2 and 3 respectively, together with the fact that 
the standard bounded extension operator i? 1 (fi) — > i? 1 (R Ar ) is also bounded in 
I? norm, to infer that for all g E H 1 ^), 

\\g\\h<C\\g\\ L 2\\g\\ Hl (N=2), (150) 

USUI. < \\9\\% 3 \\9\\% 3 < C||fllU»||s||^ (N = 3). (151) 
Using that 1 (SI) embeds into L 4 and L 6 , these inequalities lead to the estimates 

'K|| 2 4 ||V^|| 2 4 < C||u|| ia ||Vfi n ||i a ||ViZ n || Jf i (N = 2), 
Kllia||Vu"|| 2 3 < C\\Vu n \\ 3 L2 \\Vu n \\ H i (N = 3). 



|tT • Vu"| 2 < 

n 



(152) 

By the elliptic regularity estimate ||Vif||iji < ||u|| ff 2 < C||A-u||, we conclude 

(£ 2 ||Au"|| 2 +4CK 1 ||u n || 2 ||V?/ n || 4 (JV = 2), 
|K-VzT"|| 2 <i " " 2 V ' (153) 

[£ 2 ||Au"|| 2 + 4Ce 2 " 1 ||Vu™|| 6 {N = 2 or 3). 

for any £ 2 > 0. Plug this into (|147|) and take e\, £ 2 > satisfying v — e\ > 
and e := v — e\ — v(5 — 8e 2 /ei > 0. We get 

^(||ViT +1 || 2 - ||Vu"|| 2 ) + (i/-£ i )(||AiT+ 1 || 2 - ||Au"|| 2 ) +e||A?T|| 2 
< — liril 2 + — l|V?T|| 6 + uC \\ViP\\ 2 . (154) 

£l £i£ 2 

A simple discrete Gronwall-type argument leads to our main stability result: 
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Theorem 5 Let Q be a bounded domain ir 



pTV 



(N — 2 or 3) with C 3 bound- 
ary, and assume f € L 2 (0,T; L 2 (il,M. N )) for some given T > and t2° £ 
H%{n,R N )nH 2 (n,R N ). Consider the time-discrete scheme fJMj-(TJ^. Then 
there exist positive constants T* and C 3 , such that whenever nAt < T* , we have 



sup ||W fc || 2 +J2 \\&u k \\ 2 At < C 3 



E 

fc=0 



0<fc<n 



k=0 



At 



\\u k ■ Vu k \\ 2 }At<C 3 



The constants T* and C3 depend only upon Q, v and 



M a :=\\Vu \\ 2 + vAt\\Au <> \\ 2 + C \\f\\ 

Jo 



Proof: Put 

^ = ||V^|| 2 + - eOAiHA^H 2 , w n =e\\Au n \\ 2 , b n =\\f 
and note that from i|143|) we have that as long as nAt < T, 

i-T 



n 1 1 2 



1* ^ />i 

k=0 Ja 



dt. 



Then by fTSIj) . 



z n+ i + w n At < z n + CAt(b n + z n + z 3 ), 



(155) 
(156) 



(157) 



(158) 



(159) 



where we have replaced max{8/ei, 32C/ (£122), vCp} by C. Summing from to 
n — 1 and using (|158fl yields 



Zn + w k At < CM o + CAtJ2( z k + 4) =■ Vn 

k=0 k=0 

The quantities y n so defined increase with n and satisfy 

Vn+i - Vn = CAt(z n + zl) < CAt{y n + y 3 ). 



(160) 



(161) 



Now set F(y) = ln(^l + y 2 /y) so that F'(y) = -(y + y^y 1 . Then on (0,oo), 
F is positive, decreasing and convex, and we have 

F(y n+ i) - F(y n ) = F'(U)(Vn+i - Vn) > - V TV V 3 n > ( 162 ) 



whence 



Vn + Vn 



F{y n ) > F(y ) - CnAt = F{CM ) - CnAt. 



(163) 
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Choosing any T* > so that C* := F(CM ) - CT* > 0, we infer that as long 
as nAt < T* we have y n < F _1 (C*), and this together with l|160|l yields the 
stability estimate <|155fl . 

Now, using ljl53J) and elliptic regularity, we get from l|155|l that 

n n n 

^2\\u k - X7u k \\ 2 At < cY^\\Vd e \\l*\\Vu k \\w At < C J2 ll A " fe H 2Ai < c - 

k=0 k=0 k=0 

(164) 

Then the difference equation l|139|) yields 

n-l 

E 

fc=0 

This yields (|156|l and finishes the proof of the Theorem. □ 



7«+i 



At 



At < C. 



(165) 



5 Existence and uniqueness of strong solutions 

The stability estimates in Theorem|Sllead directly to the following existence and 
uniqueness theorem for strong solutions of the unconstrained formulation J^J of 
the Navier-Stokes equations. Regarding the constrained Navier-Stokes equa- 
tions there are of course many previous works; see |Am| for a recent comprehen- 
sive treatment. For unconstrained formulations of the Navier-Stokes equations 
with a variety of boundary conditions including the one considered in the present 
paper, Grubb and Solonnikov |(JSll IGS2| lay out a general existence theory in 
anisotropic Sobolev spaces using a theory of pseudodifferential initial-boundary- 
value problems developed by Grubb. 

Theorem 6 Let fl be a bounded domain in M 3 with boundary T of class C 3 , 
and let f G L 2 (0, T; L 2 (fi, M. N )), u in e H&(Cl,R N ). Then, there exists T* > 
depending only upon Q, v and M\ := \\^u- ln \\ 2 + ||/|| 2 , so that a unique strong 
solution of j i,^6)) exists on [0,T*], with 

u e L 2 (0, T*;H 2 (n, R N )) n ff^O, T*;L 2 {n, R w )), 
Vp = Vp E + vVp s E L 2 (0,T*; L 2 (n,R N )), 

where p E and p s are as in and Moreover, u £ C([0, T*], H 1 ^, M. N )), 

and V • u 6 C°°((0, T*], C°°(0)) is a classical solution of the heat equation with 
no- flux boundary conditions. The map 1 1— > || V • u\\ 2 is smooth for t > and we 
have the dissipation identity 

|i||V.«|| 2 + H|V(V.«)|| 2 = 0. (166) 

Proof of existence: We shall give a simple proof of existence based on the 
finite difference scheme considered in section using a classical compactness 
argument |Tall IT^Tl ILM| . However, in contrast to similar arguments in other 
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sources, for example by Temam |Telj for a time-discrete scheme with implicit 
differencing of pressure terms, we do not make any use of regularity theory for 
stationary Stokes systems. 

First we smooth the initial data. Given Ui n £ Hq(Q,M. ) and At > 0, 
determine u° in n H 2 (Q,R N ) by solving (I - At A) it = u in . An energy 
estimate yields 

||Vu°|| 2 + At||Au°|| 2 < ||Vui n ||||Vu°|| < ||Vu in || 2 . 

Then \\At Au°|| 2 = O(At) as At ->■ 0, so if ->■ u in strongly in L 2 and weakly 
in H . The stability constant C3 in Theorem is then uniformly bounded 
independent of At. 

We define the discretized solution u 11 by H139Jl - 1142|) of sectional and note 

jfn+l _ „yn 

— + V{u n • Vif - /" - vAu") = i/A(tT +1 - fl") + i/VV • if 1 . (167) 

With t n — nAt, we put UAt(tn) — U^(t n ) — u n for n — 0,1,2,..., and 
define UAt(t) and U& t (t) on each subinterval [t n ,t n + At) through linear inter- 
polation and as piecewise constant respectively: 

u At (t n + s) = u n + s i — j, sG[0,At), (168) 

UAt(t n + s)=u n , se[0,At). (169) 

Then (|167|) means that whenever t > with t ^ t n , 

d t UAt+V(U A t ■ VC7 Af - fAt - vAUai) = vA(U At (- + At) - U At ) + fVV • U At , 

(170) 

where f At {t) = f n for t e [t n ,t n + At). 

We will use the simplified notation X(Y) to denote a function space of the 
form X([Q,T*],Y(Cl,R N )), and we let Q = fl x [0,T*] where T* is given by 
Theorem|S| The estimates in Theorem[!)]say that UAt is bounded in the Hilbert 
space 

V := L 2 (H 2 nH^HH^L 2 ), (171) 

and also that UAt is bounded in L 2 (H 2 ), uniformly for At > 0. Moreover, 
estimate <|155|) says UAt is bounded in ClH 1 ). This is also a consequence of the 
embedding V ^ C{H r ), see [ED P- 42] or Ev, p. 288]. 

Along some subsequence Atj — > 0, then, we have that UAt converges weakly 
in Vq to some u S Vq, and UAt and UAt(- + At) converge weakly in L 2 (H 2 ) 
to some U\ and U2 respectively. Since clearly Vq <—* H 1 (Q), and since the 
embedding H 1 (Q) L 2 (Q) is compact, we have that UAt — * u strongly in 
L 2 (Q). Note that by estimate l|156|) . 

n-1 

\\uAt - ^Atlli-(Q) < II^AtO + At) - &At\\h(Q) = H^" +1 - ^H 2 ^ " CA<2 ' 

k=0 

(172) 
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Therefore £/a*( - + At) and C/a* converge to u strongly in L 2 (Q) also, so U\ = 
U 2 =u- 

We want to show u is a strong solution of 1134(1 by passing to the limit in 
1(170(1 . From the definition of /" in 1(143(1 . it is a standard result which can be 
proved by using a density argument that 

||/-/A t ||i= (o3 ^0 asAt^O. 

We are now justified in passing to the limit weakly in L 2 (Q) in all terms in 
1(17011 except the nonlinear term, which (therefore) converges weakly to some 
w G L 2 {Q). But since VE/a* converges to Vfi weakly and C/a* to u strongly 
in L 2 (Q), we can conclude C/a* ■ VE/a* converges to w • Vw in the sense of 
distributions on Q. So w = u ■ Vw, and upon taking limits in ((170(1 it follows 
that 

d t u + V(u- Vu- f- vAu) = i/VV • u. (173) 

That is, u is indeed a strong solution of ^ 134(1 . That u(0) = U[ n is a consequence 
of the continuity of the map u — > u(0) from Vb through C(i? 1 ) to i/ 1 (il,R Ar ). 
It remains to study V • u. Dot (fTT^) with V0, G H 1 ^). We get 

/ d t u-V<f> = v I V{V-u)-V4>. (174) 

This says that w = V • u is a weak solution of the heat equation with Neumann 
boundary conditions: 

dtw — vAw in Q, n-Vw = on F. (175) 

Indeed, the operator A := vA defined on L 2 (il) with domain 

D(A) = {w G H 2 {Q) | n- Vw = on T} (176) 

is self- adjoint and non-positive, so generates an analytic semigroup. For any 
4> G D(A) we have that t i— > (w(t),(f>) = — (u(i),V0) is absolutely continuous, 
and using 1(174(1 we get (d/dt){w{t), <p) = (w(t), A<j)) for a. e. i. By Ball's charac- 
terization of weak solutions of abstract evolution equations |B~a| . uj(t) = e A *w(0) 
for all i G [0,T*]. It follows w G C([0, T*], L 2 (ft)), and w(t) G D(A m ) for every 
m > Pa, theorem 6.13]. Since A m w(t) = e A ( t ~ T ^ A m w(T) if < r < t we infer 
that for < < < T*, it)(t) is analytic in t with values in D(A m ). Using interior 
estimates for elliptic equations, we find w G C°°((0, T*], C°°(Q,)) as desired. 
The dissipation identity follows by dotting with w. 
This finishes the proof of existence. □ 
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Proof of uniqueness: Suppose u\ and u 2 are both solutions of l|134|l - H136[) 

belonging to V . Put u = Ui~u 2 and Vp s = (I-V)(A- VV-)u. Then u(0) = 
and 

d t u + V(ui- Vu + u ■ Vu 2 ) = vAu ~ vVp s . (177) 
Dot with — Au and use Theorem ^ to get 

(i/Atf-i/Vftj-AtZ) < -^||Au|| 2 + ^||Vp s || 2 < -^-\\Au\\ 2 + C\\\7u\\ 2 . (178) 

Next, use the Cauchy-Schwarz inequality for the nonlinear terms, estimating 
them as follows in a manner similar to H15()|) - (|152|) . using that u\ and u 2 are a 
priori bounded in H 1 norm: 

\\ui ■ Vu\\ \\Au\\ < C||Wi||||W|| 1/2 ||Au|| 3 / 2 < e||Au|| 2 + C||Vu|| 2 , (179) 
\\u-Vu 2 \\\\Au\\ <C\\Vu\\\\Vu 2 \\ H i\\Au\\ < e\\Au\\ 2 + C\\Au 2 \\ 2 \\Vu\\ 2 . (180) 

Lastly, since u £ Vb we infer that (d t u, — Au) € i 1 (0,T) and t i— > ||Vw|| 2 is 
absolutely continuous with 

(d t u,-Au) = ~\\Vu\\ 2 . (181) 

This can be shown by using the density of smooth functions in Vb; see |Evl p. 
287] for a detailed proof of a similar result. 

Through this quite standard-style approach, we get 

^|| W|| 2 + a\\Au\\ 2 < C(l + ||AiT 2 || 2 )||Vu|| 2 (182) 

for some positive constants a and C. Because || Au2|| 2 G £ 1 (0, T), by Gronwall's 
inequality we get ||Vu|| = 0. This proves the uniqueness. □ 

Since the interval of existence [0,7^] depends only upon Mi, in standard 
fashion we may extend the unique strong solution to a maximal interval of time, 
and infer that the approximations considered above converge to this solution up 
to the maximal time. 

Corollary 2 Given the assumptions of Theorem^ system \134)) - ]lSb)) admits 
a unique strong solution u on a maximal interval [0, T max ) with the property that 
if T max < T then 

ll^Wllff 1 — v 00 ast^T niax . (183) 

For every T S [0,T max ), the approximations ma* constructed in $168\) converge 
to u weakly in 

l 2 ([o, f],H 2 n H%(n, r n )) n h\[o, f],L 2 (n, m N )) 

and strongly in L 2 ([0,f] x Q,R N ). 
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6 Unconditional stability and convergence for 
C 1 /C° finite element methods without inf-sup 
conditions 

The simplicity of the stability proof for the time-discrete scheme in section 4 
allows us to easily establish the unconditional stability and convergence (up to 
the maximal time of existence for the strong solution) of corresponding fully 
discrete finite-element methods that use C 1 elements for the velocity field and 
C° elements for pressure. To motivate the discretization, we write the uncon- 
strained Navier-Stokes formulation I jlljl in weak form as follows, in terms of 
total pressure p — p E + vp s : 

{u t + Vp- vAu + u-Vu- f,Av) =0 Vw e H 2 nH^(Q,R N ), (184) 

(Vp + vVV -u-vAu + u-Vu- f,V<p) =0 VcfreH 1 ^). (185) 

We suppose that for some sequence of positive values of h approaching zero, 
Xh C H 2 n Ho(Q,R ) is a finite-dimensional space containing the approximate 
velocity field, and suppose Yh C if 1 (il)/R is a finite-dimensional space contain- 
ing approximate pressures. We assume these spaces have the approximation 
property that 

WeH 2 nH^(n,R N )yh3v h eX hl \\A(v -v h )\\ -> as/i->0, (186) 

€ ^ity/R V7i 3cf> h e Y h , ||V(>- <t> h )\\ -> as/i->0. (187) 

As we have emphasized in the introduction to this paper, we impose no inf-sup 
condition between the spaces Xh and Yh- (We remark that in general, practical 
finite element methods usually use spaces defined on domains that approximate 
the given J7. For simplicity here we suppose f2 can be kept fixed, such that 
finite-element spaces Xh and Yh can be found as described with C 1 elements 
for velocity and C° elements for pressure. Though generally impractical, in 
principle this should be possible whenever Q has a piecewise polynomial C 3 
boundary.) 

We discretize (I184f) - H185|) in a straightforward way, implicitly only in the 
viscosity term and explicitly in the pressure and nonlinear terms. The resulting 
scheme was also derived in |JL| and is equivalent to a space discretization of 
the scheme in (|139|) - (I143|I . Given the approximate velocity wjj at the n-th time 
step, we determine p\ G Yh and € Xh by requiring 

(Vpl + i/VV ■ i% - uAul + W£ - f n , V<p h ) = \/cf> h e Y h , (188) 
Vi7™ +1 — X7iP l 

™h) + <^A< 1+1 , A^) = (VK + u n h -Vu)\ - r,Av h ) 

W h eX h . (189) 

Stability. We are to show the scheme above is unconditionally stable. First, 
we take <fih = Ph in Ijl88(l . Due to the fact that 

(p(A-vv-KWK) = o, 
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we directly deduce from the Cauchy-Schwarz inequality that 

\\Wh\\ < \\^p s «)\\ + K-vts - rii (wo) 

where 

Vp s «) = (J-P)(A-VV-K (191) 
is the Stokes pressure associated with u£. (Note Vp s (w^) need not lie in the 
space Yh). Now, taking Vh = in 1189(1 and arguing just as in l|144fl . we 

obtain an exact analog of 11451) . namely 

^(||V< +1 || 2 - ||V<|| 2 ) + {v- £l )||A< +1 || 2 

< ^ (liril 2 + IK • W£|| 2 ) + H|VpsK)|| 2 . (192) 

Proceeding now exactly as in section 4 leads to the following unconditional 
stability result. 

Theorem 7 Let £1 be a bounded domain in M. N (N = 2 or 3) with C 3 boundary, 
and suppose spaces X h C H 2 n ff^ft,!^), C ff^/M satis/y f7gg)) -f?57P. 

Assume / G £ 2 (0, T; L 2 (f2, R )) /or some given T > and S X^. Consider 
the finite- element scheme with TTien £/iere exist positive con- 

stants T* and C^, such that whenever nAt < T*, we /rave 

sup ||V^|| 2 + V ||A4|| 2 At < C 4 , (193) 

n / 7.. / „ ' ^ 

+ Il4-V^|| 2 j A<<C 4 . (194) 

Tne constants T* and C4 depend only upon f2 ; j/ and 

M , l :=||V^|| 2 + ^Ai||A ? 2O|| 2 + / T ||/T- 

Jo 

Convergence. We prove the convergence of the finite-element scheme de- 
scribed above by taking h — ► to obtain the solution of the time-discrete scheme 
studied in section, then At; — * as before. Because of the uniqueness of the so- 
lution of the time-discrete scheme and of the strong solution of the PDE, it 
suffices to prove convergence for some subsequence of any given sequence of val- 
ues of h tending toward 0. The bounds obtained in Theorem make this rather 
straightforward . 

Fix At; > 0. The bounds in Theorem[7|and in l|19(J|) imply that for all positive 
integers n < T*/At, the u% are bounded in H 2 n Hq(CI,'R n ) and the VpJJ are 
bounded in L 2 (tt,Wi N ) uniformly in h. So from any sequence of h approaching 
zero, we may extract a subsequence along which we have weak limits 

u r ^u ri in H 2 (Q,R N ), Vpl -» Vp n , u n h ■ Vu% w n in L 2 (Q,R N ) (195) 



E 

fc=0 



At 
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for all n. Then u% — > u n strongly in Hq (CI, E ) and so w n — u n ■ ViT 1 since the 
nonlinear term converges strongly in L . 

Now, for any v € i/ 2 n -ffp and € H 1 ^), by assumption there exist 

v h £ X h , 4> h e H l (9) such that if/, -> v strongly in H 2 (Cl,R N ) and V<j>h -* V<^> 
strongly in L 2 (ri,M Ar ). Applying these convergence properties in l|188|l - (|189JI 
yields that the weak limits in (|195|l satisfy 

(Vp n + z/VV-tT l -^Au n + tt n -Vu n -/",V0) = 0, (196) 

,7^+1 _ jjn 

( — vAt? l+1 + Vp" + iT-Vu" - f n , Av) = 0. (197) 

But this means exactly that u" 1 satisfies <|139[) with p n — p™ + i/p™ , where and 
are given by H14U|) - (|141I) . So in the limit /i^Owe obtain the solution of the 
time-discrete scheme studied in section Q. Then the limit At — > yields the 
unique strong solution on a maximal time interval as established in section [SJ 



7 Semigroup approach for the homogeneous lin- 
ear case 

There are many other approaches to existence theory for the Navier-Stokes 
equations, of course — Galerkin's method, mollification, semigroup theory, etc. 
We will not discuss any of them here, except to note that the linearization of the 
unconstrained system i|Tl"|) can be treated easily by analytic semigroup theory 
using Theorem ^ Take v = 1 without loss of generality, and consider Ijllfl 
without the nonlinear and forcing terms, i.e., consider the unconstrained Stokes 
equation 

u t -Aw + Vp s = (t>0, x£fl), (198) 

with the no-slip boundary condition (|135l) and initial condition (|136|) . where 
Vp s is given by (|10[) as before. In the space X = L 2 (Cl,R N ) define operators 
B and J3i by 

B u=-Au, B x u = Vp s = (I- V)Au - VV • u, (199) 

with domain D(B ) = D{B X ) = H 2 n H^(Cl,R N ). Then B is a positive self- 
adjoint operator in X with compact resolvent, and by using Theorem^] together 
with the interpolation estimate 

||Vu|| < e||At?|| + C e ||u|| 

valid for any e > for all u £ D(Bq), we deduce that 

||Biu|| < a\\B u\\ +K\\u\\ (200) 

for all u G D(Bq), where a and K are positive constants, with a < 1. 
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Theorem 8 The unconstrained Stokes operator B = Bq + B\ in the space 
X = L 2 (£l,M. N ) is sectorial and generates an analytic semigroup. The resolvent 
of B is compact and the spectrum of B consists entirely of isolated eigenvalues 
of finite multiplicity, all of which are positive. Moreover, for any a > 0, given 
M; n £ D(B a ) equation 1198$ has the solution 

u = e- Bt u in e C{[0,T},D(B a ))nC°°((0,T} : D(B m )) 

for any T > and all m > 0, and this is the unique weak solution of dtu + Bu = 
0, u(0) = Ui n in the sense of Ball JfBal/ . 

Proof. That B is sectorial is a consequence of <|200[) and the self-adjointness of 
Bq. Indeed, by a theorem on the perturbation of sectorial operators |Hel p. 19, 
theorem 1.3.2], it suffices to show that for some <f>o < 7r/2, 

a sup \\B {\- Bq^W < 1 (201) 

where So C C is the sector where 4>o < | argA| < n. By expanding any element 
of (complexified) X with respect to an orthonormal basis of eigenfunctions of 
B , for any A ^ cr{B ) we get 

\\B {\-B )- 1 \\= sup 

Fix a E (a, 1). For any /i > 0, we have \fi\ < |A — /i| whenever 3?A < 0, and 
it is straightforward to check that whenever -JJA > and |3A| > 5|A|, then 
2|/i| < I A — /-tj . Then (|2()l(l follows, proving that B is sectorial. 

That (A— B)^ 1 is compact for A ^ o-(B)Ua(Bo) follows from the compactness 
of (A — Bq)^ 1 together with the identity 

(A - B)- 1 = (A - Bo)- 1 + (A - B )- 1 B 1 (X - B)~ l . 

It follows that the spectrum of B is discrete, consisting only of isolated eigen- 
values of finite multiplicity |Kal III. 6. 29]. 

Suppose now that (A — B)u = for some non-zero u € D(B), so Xil = 
—VA.H — VV • u. Then the function w — V ■ u satisfies Xw = —Aw in fi, 
ft ■ Vw = on T, i.e., (A + A)w = (see ljT7S))). So if A ^ o-(-A) C K+, then 
V • u = 0, and since ft ■ u = on T we have u = Vu. Then 

X(u,u) = ( - VAu,u) = ( - Au,Vu) = ||Vu|| 2 , 

so A > 0. If A = 0, then V • u is constant, but J n V- , u = 0soV-'U = and 
arguing as above we infer u = 0. Hence is not an eigenvalue, and so is in 
the resolvent set of B. 

Lastly, for any a > 0, given u m € D(B a ), the regularity results for e~ Bt u ln 
are standard consequences of the fact that B a is an isomorphism between its 
domain and X and commutes with e~ Bt |Pal p. 74, Theorem 6.13]. For unique- 
ness, see |Baj . □ 



A — /i 
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Remark 3. The equation Bit — f has an interesting interpretation in terms of a 
stationary Stokes system with prescribed divergence. Given any / £ L 2 (fl,R. N ) 
there is a unique u £ H 2 n Hq(Q,M. n ) such that Bu = /, since is in the 
resolvent set of B by the above theorem. We can write Vf = f + Vq where 
q £ H 1 ^) with / n q = 0. Since Bu = -VAu - VV • u, we have Vf = -VAu, 
so V(g + V • u) — 0. Let p s be the Stokes pressure associated with u. Then 
(u,p s ) form a solution to the Stokes system 

-Au + Vp s = f in ft, (202) 
-V ■ u = q in n, (203) 
u = on T. (204) 

As a corollary, we can characterize the domains of positive integer powers of 
B by using the regularity theory for the stationary Stokes equation (see for 
example Soh p. 123, theorem 1.5.3] or [Tell p. 23, proposition 2.2]). 

Corollary 3 Let Q, be a bounded domain with C 2m boundary Y , where m > 1 
is an integer. Then 

D{B m ) = {u | u £ H 2m (Cl,R N ), u = Bu= ... = B m - 1 u = onT}. 

Proof: When m — 1, the conclusion is true. Suppose it is true when m = 
k — 1. When m = k, take any u £ D(B k ). By the definition of D(B k ), we 
have u £ D(B k ^ 1 ) and Bu £ D(B k ^ 1 ). By assumption, / := Bu £ H 2k - 2 
and B k ~ lr u = B k - 2 (Bu) = on T. Since V is bounded on H 2k ~ 2 JTHTJ I, 
Remark 1.6] we find that q £ H 2k - 1 (Q). Now <j2"U2"j) - (l2"U4l) hold, and we can use 
the regularity theory of the stationary Stokes equation cited above to conclude 
u £ H 2k (fl,R N ). This finishes the proof. □ 

1 /2 

Remark 4. We note that B and Bo have the same domain and that D(B ) 
is the closure of D(Bq) — H 2 n Hq(£1, M. n ) in norm equivalent to 

\\u\\ 2 xl/2 = \\Bl /2 u\\ 2 = ( - Au,u) = ||W|| 2 , 

the ordinary H 1 norm. So D(Bl /2 ) = H^(ft,R N ). It is known that if B 
has bounded imaginary powers then for < a < 1, D(B a ) can be obtained 
by interpolation between X and D(B) = D(B ) and so D(B a ) = D(B%). 
The result that indeed B + cl has bounded imaginary powers for some c > 
apparently follows from a recent analysis of Abels |Ab| related to the formulation 
of Grubb and Solonnikov (although the final result in |Ab| is stated in terms of 
the constrained Stokes operator in divergence-free spaces). 

8 Non- homogeneous side conditions 

Looking back at the Stokes pressure p s associated with it, one recognizes that the 
no-slip boundary condition for u was essential for getting the crucial equalities 
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l|$§l) - (|lUlfl using Lemma|3] So the important question arises, if general boundary 
conditions u = g on T are imposed, do we still have an unconstrained formulation 
like l|134|) - (|136|) ? Moreover, what can we say if the velocity field is not divergence 
free but is specified as V • u = hi Such issues are likely to be relevant in the 
analysis of problems involving complex fluids and low Mach number flows, for 
example. 

In this section we develop and study an unconstrained formulation for such 
non- homogeneous problems. In this new formulation, V -u — h satisfies the heat 
equation with no-flux boundary conditions. The main theorem of this section 
establishes existence and uniqueness for strong solutions. 



8.1 An unconstrained formulation 

Consider the Navier-Stokes equations with non-homogeneous boundary condi- 
tions and divergence constraint: 

d t u + u-Vu + Vp = vAu + f (t>0,x€Cl), (205) 

V-u = h (t>Q,xefl), (206) 

u = g {t>Q,xeT), (207) 

u = u in (i = 0,x6Sl). (208) 

What we have done before can be viewed as replacing the divergence constraint 
l|2Utjl) by decomposing the pressure via the formulae in © and lfTU)l in such a 
way that the divergence constraint is enforced automatically. It turns out that 
in the non-homogeneous case a very similar procedure works. One can simply 
use the Helmholtz decomposition to identify Euler and Stokes pressure terms 
exactly as before via the formulae © and (|1C)|> , but in addition another term is 
needed in the total pressure to deal with the inhomogeneities. Equation 10 is 
replaced by 

d t u + V(u-Vu - /- uAu) + Vp gh = iA7(V • u). (209) 

The equation that determines the inhomogeneous pressure p g h can be found by 
dotting with for (ft e H 1 ^), formally integrating by parts and plugging in 
the side conditions: We require 

(Vp fffc ,V0) = -(d t (ft- g),<ft) r + {d t h,4>) + (uVh,V4>) (210) 

for all (ft e H 1 ^). With this definition, we see from l|2U9|l that 

(d t u, V(ft) - (d t (rt ■ g), 4>) r + (d t h, (ft) = -u- h),V(ft) (211) 

for every (ft S H 1 ^). This will mean w := V • u — h is a weak solution of 

d t w = vAw in Q, n ■ Vw = on T, (212) 

with initial condition w = V • u- m — h\t—o 4 ^° ^ ne divergence constraint will be 
enforced through exponential diffusive decay as before (see (|232ll below). 
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The total pressure in l|2U5l) now has the representation 

p = p E + VPs + Pgh, (213) 

where the Euler pressure p E and the Stokes pressure p s are determined exactly 
by JHJ and l|l Ojl as before, and p g h is determined up to a constant by the forcing 
functions g and h through the weak- form pressure Poisson equation H21()[l . (See 
Lemma El below.) Our unconstrained formulation of i|2U5[) - i|2U8|) then takes the 
form 

d t u + u-X7u + Vp E + v\7p s + Vpgh = vAu + f (t > 0, x S fl), (214) 

u = g (t > 0, x e P), (215) 
u = u- m (t = 0, x G fi). (216) 

Although the definition of Stokes pressure does not require a no-slip velocity 
field, clearly the analysis that we performed in section 2 does rely in crucial ways 
on no-slip boundary conditions. So in order to analyze the new unconstrained 
formulation, we will decompose the velocity field u in two parts. We introduce 
a fixed field u in x [0, T] that satisfies u = g on T, and let 

v = ti-u. (217) 

Then v = on T. With this v, similar to © and (|ll)f> we introduce 

Vq E = (V -I)(v-Vv- f), Vq s = (I-V)Av- VV • v. (218) 

Then we can rewrite ()214ll as an equation for v: 

d t v + v-Vv + Vq E + vVq s + V{u ■ Vv + v ■ Vu) = vAtf+f- f, (219) 

where 

/ := d t u + V(u • Vw - vAu) - vW ■ u + Vp gh . (220) 

8.2 Existence, uniqueness and dissipation identity 

We will first answer questions concerning the existence and regularity of u and 
Pgh, then state an existence and uniqueness result for strong solutions of the 
unconstrained formulation (|214|l - (|216|l . Let be a bounded, connected domain 



in M. N (N — 2 or 3) with boundary T of class C 3 . We assume 

^^(O,!^, (221) 

/ e H f := L 2 {0,T;L 2 (n,R N )), (222) 
ge H g := H 3/4 (0,T; L 2 {T,R N )) n L 2 (0,T; H 3/2 (T,R N )) 

n {g | d t (n ■ g) e i 2 (0, T; #- 1/2 (r))}, (223) 

ftG-ffft :=i 2 (0,T;iJ 1 (fi))ni7 1 (0,T;(ff 1 )'(fi)). (224) 
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Here (H 1 )' is the space dual to H 1 . We also make the compatibility assumptions 

g = M in when t = 0, x e T, (225) 
<a t (n-3),l) r = (5^,l) . (226) 

We remark that most of the literature on nonhomogeneous Navier-Stokes prob- 
lems |Lal ISoll IGrll IOS1I IOS2j treats the constrained case with h = in SI and 
imposes the condition n ■ g — on T. Amann recently studied very weak so- 
lutions without imposing the latter condition, but only in spaces of very low 
regularity that exclude the case considered here Am2 . 
We define 

V := L 2 (0,T;H 2 (Q,R N )) n H 1 (0,T; L 2 (Sl,R N )), (227) 

and note we have the embeddings f [TaTl p. 42], [Ev| p. 288], |TeTl p. 176]) 

V^C^TIH 1 ^,^)), H h ^C([0,T},L 2 (n)). (228) 

Notice that we have always used an arrow or tilde to denote a vector. So, without 
confusion, we can use Y(fl) to denote Y(f2, R-^) or Y(p,) as appropriate, and 
further use X(Y(9,)) to denote X(0,T;Y(Q)). 

Lemma 6 Assume \221\) - \22b)) . Then, there exists some u £V that satisfies 

u(0)=tt- m , u\ r =g, (229) 
and there exists p g h £ L 2 (H 1 (fl) /M.) satisfying ]210\) . Moreover, 

INly ^ ^(ll5llff3/4(L 2 (r))nff 3 /2(r)) + Il"in||ffi(f2)), (230) 
\\Pgh\\L^(m(n)/M) < C(\\d t (n ■ g)\\L*(H-i/2(r)) + II^IU 2 (ff 1 )nffi((,F/i)'))- ( 231 ) 

Proof: (i) By a trace theorem of Lions and Magenes |LM1 vol II, Theorem 2.3], 
the fact g € H 3 ^{L 2 {T)) n L 2 (H 3 / 2 (Y)) together with &ZH and the compati- 
bility condition (|225|l implies the existence of u <E V satisfying Ij229|l . 

(ii) One applies the Lax-Milgram lemma for a.e. t to (|210|l in the space of 
functions in 1 (17) with zero average. We omit the standard details. □ 

Theorem 9 Let Q be a bounded, connected domain in M. N (N = 2 or 3) and 
assume \221}) - $226]) . Then there exists T* > so that a unique strong solution 
of ggp- fH^) exists on [0,T*], with 

u GL 2 (0, T*;H 2 (Q,,R N )) n ff x (0, T*;L 2 (fl, R N )), 
p =vp s + p E + Pgh e L 2 (0, T*;H 1 (Q)/R), 

where p E and p s are defined in and after introducing the u and p g h from 
Lemma\B Moreover, u <E C([Q, T*] , H 1 (fi, R N )) and 

V -u-he L 2 {0, T*; H^Sl)) H H\0, T*; (H 1 )'^)) 
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is a smooth solution of the heat equation for t > with no-flux boundary condi- 
tions. The map 1 1— * || V • u — h\\ 2 is smooth for t > and we have the dissipation 
identity 

-7™||V -it - /i|| 2 + H|V(V •w- /i)|| 2 = 0. (232) 

If we further assume h £ H h . s := L 2 (0, T; H 2 (fl))r\H 1 (Q, T; L 2 (fl)) andV -u in £ 
H 1 ^), then 

V ■ u £ L 2 {0, T*;H 2 (Q)) n £^(0, T*;L 2 (Vt)). 



Proof: First rewrite (|214|l as (|219|l . Then we note that there are only two 
differences between (I219|) and (|134l) : 

(i) There is an extra forcing term / in i|219l) . But by Lemma HJ all terms in 
/ are known to be in L 2 (L 2 (Vt)) and thus they won't be a problem. 

(ii) Equation (|219|l has some extra linear terms: 

V{u-\7v + v-\7u). (233) 

We know u G V C([0, T], i? 1 (^, 1^)), so we can discretize these terms 
explicitly by setting u n — u(nAt). Similar to l|153[l . we get 

llP^.V^II^ellA^p + lllSll^llV^II 2 . (234) 

We estimate the other term in Ij233(l by using Gagliardo-Nirenberg inequalities 
|Frl Thm. 10.1] and the Sobolev embeddings of H 1 into L 3 and L 6 : 

Jq|A,i|^ /2 ||,i|^ 2 <q|A,i|V 2 ||v,i|V 2 {N = 2), 

\\v < \ (235) 

^ciiAirii 1 / 2 !!^!^/ 2 < ciiAvii 1 / 2 !^^!! 1 / 2 (jv = 3). 

Then for N = 2 and 3 we have 

\\r(v^u)f < N| 2 l~I|V«|| 2 < e\\Av\\ 2 + lll^HWH 2 . (236) 

With these estimates, the rest of the proof of existence and uniqueness is essen- 
tially the same as that of Theorem^ and therefore we omit the details. 

To prove the regularity of V • u, we argue in a manner similar to the proof 
of Theorem We go from (|214|l to (|209l) by using |JSJ and (fTOft . Then using 
(j2TU|) we get (f^TTf) for any cj> £ H 1 ^). With w = V -u-h, taking (f> £ D(A) as 
in (|176|) . we have 

(w, 0} = (n • g, 4>) r - (u, V<£) - (h, (f>), (237) 

therefore t i— > (.w,(p) is absolutely continuous, and (|211|l yields (d/dt)(w,cf)) — 
(w, Acf) for a.e. t. This means u; is a weak solution in the sense of Ball Ba , 
and the rest of the proof goes as before. 
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If we further assume h € H h s and V • u m € H 1 ^), then w(0) G 7J 1 (SI). We 
claim 

fl^fi) = D^-A) 1 / 2 ). (238) 
Then semigroup theory yields w € C([0, T*}, D((— A) 1 / 2 )), so since 

= ( - Aw, d t w - uAw) = ^ll Vw ll 2 + HI Aw ll 2 ( 239 ) 

for i > 0, we deduce w € L 2 (0,T*; H 2 (Q)) n F x (0, T*; L 2 (S1)), and V • u is in 
the same space. 

To prove note X := D((-A) 1 / 2 ) is the closure of D(A) from (jT75|) in 

the norm given by 

\M\ 2 x = \M\ 2 + \\{-A) l/2 w\\ 2 = {{I -vk)w,w) = f \w\ 2 + u\Vw\ 2 . 

Clearly X C i/ 1 (il). For the other direction, let w <E 1 (SI) be arbitrary. 
We may suppose w £ C°°(S7) since this space is dense in iJ 1 (SI) . Now we 
only need to construct a sequence of C 2 functions w n — ► in H 1 norm with 
n ■ Vro n = n - \7w on T. This is easily accomplished using functions of the form 
ib„(i) = ^„(dist(a;, r))n •Vw(x), where £ n (s) = £,(ns)/n with £ smooth and 
satisfying £(0) = 0, £'(0) = 1 and f(s) = for s > 1. This proves JH^. 
We can prove the uniqueness by the same method as in Theorem □ 



9 Isomorphism theorems for non-homogeneous 
Stokes systems 

Drop the nonlinear term and consider the non-homogeneous Stokes system: 

d t u + Vp-vAu = f (t>0,x£Q), (240) 
V -u = h {t > 0, x G SI), (241) 

u = g (t > o, x e r), (242) 

u = Win (i = 0, a; G SI). (243) 
The unconstrained formulation is 

d t u + Vp-vAu = f {t>0,xeQ), (244) 

u = g (t > o, x e r), (245) 

u = u- m (t = 0, x e SI), (246) 

with 

Vp = (I - V)f + vV Ps + Vp 3/1 , (247) 

where p s and are defined as before via (|10|l and (|210|l . 

The aim of this section is to obtain an isomorphism between the space of 
solutions and the space of data {/, g, h, Ui n }, for this unconstrained formulation 
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and for the original Stokes system. In examining this question we are motivated 
by the classic works of Lions and Magenes |LM| which provide a satisfactory 
description of the correspondence between solutions and data for elliptic bound- 
ary value problems. In the spirit of these results, a satisfactory theory of a given 
system of partial differential equations should describe exactly how, in the space 
of all functions involved, the manifold of solutions can be parametrized. Yet we 
are not aware of any such complete treatment of the non-homogeneous Stokes 
system. (See further remarks on this issue below.) 

First we consider the mapping from data to solution. Thanks to the absence 
of the nonlinear term, we can repeat much easier what we did in the proof 
of Theorems and and get the global existence and uniqueness of a strong 
solution of 1244(1 - (|247|l under the same assumptions as Theorem The data 
{/, g, h, Win} he inside the space 

IIf := H f x H g x H h x H um (248) 

from H221fl - 1224fl . an d need to satisfy the compatibility conditions (|225|1 - 1226() . 
Corresponding to such data, we get a unique solution u of (|244(1 - (|247|) in the 

space 

H u := L 2 (0, T; H 2 (Q, R N )) n H\0, T; L 2 (tt, R N )) 

n{u\d t (n-u)\ r e L^Q^T-H- 1 ' 2 ^))}. (249) 

The total pressure p lies in 

H p := L^Tjff^VM), (250) 

and the pair {u,p} satisfies (|240l) . (|242[) and (|243[1 . As in Theorem we can 
show w = V • m — h satisfies a heat equation with no- flux boundary conditions. 
Equation (|24f \ says that w = 0, and this will hold if and only if w(0) — 0, i.e., 
the following additional compatibility condition holds: 

V-u in = h(0). (251) 

For the non-homogeneous Stokes system (|240f) - (|243|l . then, we define the data 
and solution spaces by 

n F . c := j{/,g>,w in } G n F : (ESI and (J2SU hold j, (252) 

Uu := H u x H p . (253) 

From what we have said so far, we get a map IIi? c — > Hjj by solving the 
unconstrained system 1)24411 -- (|247|l . Due to the absence of nonlinear terms, the 
estimates in the proof ensure that this map is bounded. In the other direction, 
given {u,p} € Hu, we simply define {/, g, h,u in } using (|240|I - H243() and check 
that this lies in Hf.c- 



Divorcing pressure from viscosity 



41 



Note that in Theorem [§| one has more regularity on V • u if one assumes 
more on V • Ui n and h. Correspondingly, like Hh.s defined in Theorem |5J we 
introduce spaces of stronger regularity by 

H lan . s := H 1 (fi, R N ) n {u in | V ■ u in G H x (Sl)}, (254) 
U F . S := H f x H g x H h . s x H uin . s . (255) 

The solution u then lies in 

H u . s := L 2 (0,T;H 2 (n,R N )) n tf^O, T; L 2 (ft, M^)) 

n{u| V-tZGi 2 (0,T;i? 2 (fi))niJ 1 (0,T; J L 2 (f7))}. (256) 

(Note, if u S iT„. s then d t u G L 2 (iJ(div; Q)) so n ■ <9 t u G L 2 {H~^ 2 {T)).) So 
as an alternative to the spaces in H252ll - (|253|) . we also obtain an isomorphism 
between the data and solution spaces with stronger regularity defined by 

n^.c.s := {{/,<?, /i,u in } G n F . s : (EH, and (gFTJ hold |, (257) 

Ilu. s ■= Hu.s x H p . (258) 

Summarizing, we have proved the following isomorphism theorem for the 
non- homogeneous Stokes system (|240[I - H243[) . 

Theorem 10 Let VI be a bounded, connected domain in M. N with N any positive 
integer > 2, and letT > 0. The map {/ , g, h,u- ln \ i— > {u, p}, given by solving the 
unconstrained system f£^| )- JI^7fl , defines an isomorphism from Hf.c onto Tljj- 
The same solution procedure defines an isomorphism from II i^. c . s onto Hu.s- 

Remark 5. For the standard Stokes system with zero-divergence constraints 
V • win = and h = 0, existence and uniqueness results together with the 
estimates 

SUP \\u(t)\\ H l + ||m|| L 2 (0!T . H 2 ) + |b|| £2(0^.^1 /R ) 

0<t<T 

< C(||/|| L 2 (0!T . L 2) + \\u iu \\ H l + ||.9|| /f 3/4 (L 2 (r)) + ||.9|| L 2 (ff 3/2 (r)) ) (259) 

were obtained in the classic work of Solonnikov [Soil Theorem 15], where more 
general L p estimates were also proved. (Also see |GS1IIGS2] .1 However, instead 
of the necessary compatibility condition 

f H -$=°> ( 26 °) 

Solonnikov made the stronger constraining assumption that both the data g and 
solution u have zero normal component on T, and correspondingly his estimates 
do not contain a term ||9t(r? • gJW^fH- 1 / 2 ^)) on the right hand side of Ij259|l . 
(Note that when V • Ui n = and h = 0, we have J r ft ■ g\t=o — Jq V • Ui n = by 
l{2"2"5|) . whence (pp is equivalent to 
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Remark 6. For the unconstrained Stokes system l|244[l - l|246|) there is an extra 
subtlety in determining an isomorphism from data to solution. We obtain a 
unique solution pair {u,p} £ Hjj given any data {/, g, h, u; n } £ lip that satisfy 
only the compatibility conditions (|225|l and (|226J) without (|251J) . Consequently 
the map from data to {u, p} is not one-to-one. And, in the other direction, given 
{u,p}, we can recover 

f = d t u + \7p-vAu, g = u\ r , u in = u\ t=0 . (261) 

But how are we to recover hi We need to use the fact, that follows from the 
definition of p g h in (|210|l . that V • u — h satisfies a heat equation with no- flux 
boundary conditions. In fact, to be able to recover h we need to know one more 
item, hi n , the initial value of h. We have 

h = V-u-w (262) 

where w is the solution of 

dtiv = vAw in Q, n ■ Vro = on T, w(0) = V • u\t=o ~ h- m . (263) 

This procedure indicates that we should count the triple {u, p, hi n } as our solu- 
tion in order to build an isomorphism with the data. Of course, the regularity 
of /iin must match that of h, recalling the embeddings in J22HI- 

Consequently, we see that solving the unconstrained system 1244|) (|247[) de- 
fines an isomorphism between the data spaces 

ftf.c := Uf,9,h,a in } £ U F : QHty and ^FM hold \, (264) 

ILf. c , s := {{/, g,h,ti in } en F , s : {229 and {22$ hold }, ( 265 ) 

and, respectively, the solution spaces for {u, p, h[ n } given by 

n UM = H u x H p x H hin , H hm =L 2 (n), (266) 

I4;7.s = H u , s X H p X Hhin.s, Hhin.s = H 1 (fi) . (267) 
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Sciences at the National University of Singapore. RLP acknowledges support 
by the Distinguished Ordway Visitors Program of the School of Mathematics, 
and the Institute for Mathematics and its Applications, at the University of 
Minnesota. 
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